NASA TECHNICAL NOTE 


NASA TN D-7834 


CO 

oo 







AMASA-TN-D-7834) APPLICATION OF A 
PARAHETEB IDENTIFICATION TECHNIQUE TO A 
HINGELESS HELICOPTER ROTOR (NASA) 68 p 
HC S4.25 CSCL 0 1 B 

V 


Bl/02 


1175 - 12906 ^ 


Unclas 

06127 



APPLICATION OF 

A PARAMETER IDENTIFICATION TECHNIQUE 
TO A HINGELESS HELICOPTER ROTOR 


by Gerd Kanning and James C. Biggers 

Ames Research Center 

field, Calif. 94035 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, 0. C. • DECEMBER 1974 



"at onal technical 
formation service 

#i %sitriww“ 









1. Report No. 

MSA TN D-7834 


2. Government Accession No, 


4. Title and Subtitle 


APPLICATION OF A PARAMETER IDENTIFICATION TECHNIQUE TO A 
HINGELESS HELICOPTER ROTOR 


7. Author{$! 

Gerd Fanning and James C. Bigger s 


9. Performing Organization Name and Address 

NASA Ames Research Center 
Moffett Field, Calif. 9k 035 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 205^6 


3, Recipient's Catalog No. 


5. Report Date 
i December 197 ^ 


6. Performing Organization Code 


8. Performing Organization Report No. 

A-5289 


10, Work Unit No. 
760-63-01 


11. Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Note 


14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 


A mathematical model of a gyro-controlled, three-bla&ed hingeless helicopter rotor was 
developed and parameters of the model were, estimated, using, a parameter identification technique. 

The flapping and feathering degrees of freedom of the blades were modeled. The equations of the 
model contain time-varying, periodic coefficients due to the forward speed of the rotor. A digital 
simulation of the analytical model was compared with wind-tunnel measurements to establish the 
validity of the model. Comparisons of steady-state and transient solutions of. the analytical model 
with the tunnel measurements gave reasonably good matching of gyro angle but less satisfactory- 
matching of hub moment measurements. Further improvements were obtained by use of a parameter 
identification technique to adjust as many as 10 parameters of the analytical model. The sensitiv- 
ity of the blade response to small changes in the parameters was also calculated. The most sensi- 
tive parameters were found to be collective- pitch , blade twist, tip loss factor, blade to gyro 
coupling constant, Lock number, and equivalent flapping hinge spring constant. 


17. Key Words (Suggested by Author(sf) 

Helicopter parameter identification 
Helicopter rotor modeling 
Mathematical modeling 
’ Helicopters 
Modeling 
Identification 

18. Distribution Statement 

Unclassified - Unlimited 

Cat. ce 

19. Security Classif . (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

) 

Unclassified 

Unclassified 



68 



For sale by the National Technical Information Service, Springfield, Virginia 22151 








TABLE OF CONTENTS 


Page 


SYMBOLS v 

SUMMARY 1 

INTRODUCTION 1 

DEVELOPMENT OF THE ANALYTICAL MODEL . 2 

Blade Equations . . , . . 3 

Gyro Equations 5 

Rotor Equations 7 

Method of Solution and Transformation of Equations to 

Measurement Coordinates , ....... * 8 

BEHAVIOR OF SIMULATION MODEL * 9 

Parameters of Model 9 

Time Response . 10 

APPLICATION OF PARAMETER IDENTIFICATION TECHNIQUE lb 

COMPARISON OF MEASURED AND SIMULATED DATA l6 

Discussion of Measured Data l6 

Steady-State Case: Locked Gyro . , . l8 

Transient Case: Free Gyro . . 20 

Discussion of Results 2h 

CONCLUSIONS 28 

APPENDIX A - DERIVATION OF BLADE EQUATION 29 

Inertia Moments * * 29 

Aerodynamic Moment . 31 

Mechanical Moments 36 

APPENDIX B - DERIVATION OF GYRO EQUATION b0 

APPENDIX C - HUB MOMENT CALCULATIONS FROM MEASURED DATA U8 

REFERENCES 51 

TABLES 53 





B^ 


$0^ 




iii 



SYMBOLS 


A(t) 

a 

B 


b 

C 


C 

C 2 


ci(^),c2(»; 
C3(-<0;c'iK*), , 
C5(^),C6(^), 
C7W ' 


time- varying matrix defined in appendix D (table l) 
lift-curve slope equal to 6.28 

frame formed by vectors (bj,b 2 »b 3 ) which defines the blade 
frame 

tip loss factor 
number of blades 

orthogonal unit vectors defined as the B frame 

frame formed by vectors (£ 1 , 02 , 03 ) which defines the gyro cage 
frame ’ ' ' 

transformation matrix between state variables and measurement 

coupling constant between blade and angular' displacement of 
gyro 

drag coefficient, — r 

(l/2)p* rrR^ 

L 

lift coefficient, — - 

(1/2)p? 2 tfR^ 

defined by equation (A35) 


n 


(£1,02,03) 


damping constant for nonrotating part of ' swashplate , 
Nm/r ad/sec (ft-lb /rad/sec ) 


thrust coefficient. 


itE 2 p(£JR b ) 2 


blade chord, m (ft) 
defined by equation (A20) 
defined by equation (A20) 

orthogonal unit vectors defined as the C frame 

damping constant about feathering axis (rotating axis), 
Nm/rad/sec (ft-lb /rad/sec) 


PRECEDING 


page blank not FILM® 


V 



(». 

e 

F *x 

F Sy 

FWD 

f(t) 

G 

g rpm 

( 5 l.g 2 .i 3 ) 


D 

f G 

X GE 

13 


I ,1 ,1 
x* y * z 

K 


K1,B2,B3, 

B4,B5 


R 


drag, N (lb) 

time rate of change with respect to inertial space 
hinge offset, m (ft) 

force at hub in rotating coordinates along 
force at hub in rotating coordinates along 
forward direction 
column vector of forcing functions for rotor equations 
frame formed by vectors (ii»g2»g3) ’which defines the gyro frame 
rotational speed of gyro, rpm 

orthogonal unit vectors defined as the G frame 

angular momentum of gyro, Nm 2 /sec ( slug-ft 2 /sec ) 

distance from hub to load cell moment sensor, m (ft); 

E LC = (°»°* h LC z ) 

distance from hub to location of shaft moment sensor, m (ft); 

5 S = (°» 0 » h s z ) 

diametric moment of inertia of gyro, kg-m 2 (slug-ft 2 ) 
gyro inertia matrix defined by equation ( B 13 ) , kg-m 2 (slug-ft 2 ) 
equivalent inertia of gyro + 1.5C2l x , kg-m 2 (slug-ft 2 ) 
polar moment of inertia of gyro, kg-m 2 (slug-ft 2 ) 
inertia characteristics of blade, kg-m 2 (slug-ft 2 ) 
stationary spring constant for swashplate, Nm/rad (ft-lb/rad) 

defined by equation (A 34 ) 

ratio of gyro angular velocity to rotor angular velocity 
spring constant of blade, Nm/rad (ft-lb/rad); 

k g = I y fl 2 [p 2 - 1 - - $ tan 63] 


X axis, N (lb) 
Y axis, N (lb) 


vi 



L 

M 

exc 

M exc xc 

^exo, 


'yc 


% 

m lc. 


M 


M, 


M 

X 

M 

y 

*v 


% 


aero 

inertia 


^inertia 

^mech 

^meeh 

M x. 


spring constant about feathering axis, Nm/rad (ft- lb/rad) 
lift, N (lb) 

excitation of the gyro by the blade in the R 
Nm (ft-lb) 

x-eomponent of M in G frame, Nm (ft-lb) 

y-component of M exc in C frame, Nm (ft-lb) 

applied moment to gyro along X in the C frame, Nm (ft-lb) 

applied moment to gyro along Y in the C frame (Nm (ft-lb) 

load cell roll moment in stationary coordinates measured 0 . 91 *+ m 
(3 ft) below hub, Nm (ft-lb) 

load cell pitch moment in stationary coordinates measured 0.91*+ m 
(3 ft) below hub, Nm (ft-lb) 

moment about X axis in the B frame 

moment about Y axis in the B frame 

aerodynamic moment about flapping axis, Nm (ft-lb) 

inertia moment about feathering axis, Nm (ft-lb) 

inertia moment about flapping axis, Nm (ft-lb) 

mechanical moment about feathering axis, Nm (ft-lb) 

mechanical moment about flapping axis, Nm (ft-lb) 

hub moment in stationary coordinates (left roll positive), Nm 
(ft-lb) 


frame along X r , 


% s 


hub moment in stationary coordinates (nose-up positive), Nm 
(ft-lb) 


M *r 

>Vr 



shaft bending moment about X r axis 0.369 m (0.917 ft) below 
hub, Nm (ft-lb) 

shaft bending moment about Y r axis 0.369 m (0.917 ft) below 
hub, Nm (ft-lb) 

roll moment at hub in rotating coordinates (positive left roll), 
Nm (ft-lb) 

pitch moment at hub in rotating coordinates (positive nose-up), 
Nm (ft-lb) 


vii 



n nth blade 

p ratio of flapping natural frequency to rotor frequency — 

q dynamic pressure, , N/m 2 (lb/ft 2 ) 

r frame formed by vectors (ri,r2,?3) which defines the shaft axis 

system rotating at rotor speed 

R. blade radius, m (ft) 

b 

r distance along blade from center of rotation, m (ft) 

(r ls r 2 ,r3) orthogonal unit vectors defined as the R frame 

S frame formed by vectors (si,s 29 S3) which defines the stationary 

axis system of the rotor; also referred to as the inertial 
frame 


(si *82,83) orthogonal unit vectors defined as the S frame 
T rotor thrust, DJ (lb) 


t time, sec 

U ,U ,U radial, tangential, and perpendicular wind velocity components, 

E T P m/sec (ft/sec) 


u input vector 

V airspeed, m/sec (ft/sec) 


X,Y,Z 


x 

X 


right-handed coordinate system with the positive X axis aft 
over helicopter tail section and the positive Z axis 
pointing up 

ratio of radial distance to blade radius, 
state variable 




y 

a = 9-4 


os 

s 

6 


ft 


P 


measurement vector 

angle between relative wind and blade chord axis - angle of 
attack, deg 

angle of attack of hub plane 

flapping angle of blade (positive up), deg 

precone angle, deg 



&1 


Y 

5 3 

«1 

S 2 

c 

8 

e 

o 

§1 


eye 


sp 


e l » e 2» e 3 
8 


1C 


8 


IS 


v 

»i 

v 2 


sp 


flapping angles of blades 1, 2, and 3, 

.4 


Lock number, y - 


:paR b 


deg 


pitch flap coupling, deg 

gyro angle, left roll positive, rotating coordinates, deg 

gyro angle, nose-up positive, rotating coordinates, deg 

lag angle of blade, deg 

blade sweep'' angle , deg 

feathering angle of blade, deg 

collective pitch, deg 

blade twist, deg 

cyclic pitch angle, deg 

swashplate pitch angle, positive nose-up (stationary coordinates), 
deg 

feathering angles of blades 1, 2, and 3, deg 

cyclic pitch, nose up at ip » 0, deg 

cyclic pitch, nose up at \p = 9C°> deg 

inflow positive up, X + X x cos ip + Q x sin ip 

ox y 

V cos a 

advan c e r at i o , 

swashplate roll moment input (stationary coordinates, positive 
left roll), Nm (ft- lb) 

swashplate pitch moment input (stationary coordinates, positive 
nose-up), B in. (ft-lb) 

air density, kg/m 3 (slug/ft 3 ) 

inflow angle, deg 

swashplate roll angle, positive left roll, stationary coordinates, 
deg 

azimuth angle measured in direction of rotation, deg 


lx 



* 


o 


£2 


^C/R 

^C/S 

“g/C 

W R/S 


pitch link angle, deg (see fig. 3) 
angular velocity of rotor, rad/sec 
blade flapping natural frequency, rad/sec 
gyro angular velocity, rad/sec 

angular velocities of blade defined by equation (A3), rad/sec 

angular velocity of C frame with respect to R frame, rad/sec 

angular velocity of C frame with respect to inertial space, 

rad/ sec 

angular velocity of G frame with respect to C frame, rad/sec 

angular velocity of R frame with respect to inertial space, 

rad/sec 


x 



APPLICATION OF A PARAMETER IDENTIFICATION TECHNIQUE 


TO A HINGELESS HELICOPTER ROTOR 
Gerd Kanning and James C. Biggers 
Ames Research Center 


SUMMARY 


A mathematical model of a three-bladed hingeless rotor 'with a high-speed 
gyro in the feathering system was developed and parameters of the model were 
estimated using a parameter identification technique. The flapping degree of 
freedom of each "blade was modeled by assuming the existence of a virtual 
hinge near the hub. The feathering degree of freedom was modeled in conjunc- 
tion with the motion of the gyro. The gyro prescribed the feathering 
sequence to each of the blades. In the aerodynamic moment terms, the reverse 
velocity effect of the retreating blades was included. The complete gyro and 
rotor equations contain time-varying periodic coefficients due to the forward 
speed of the rotor. A digital simulation of these equations was obtained to 
generate steady-state and transient solutions of the mathematical model. For 
steady-state operation, with the gyro inoperative, the model response was 
matched quite accurately to the measured data by adjusting selected model 
parameters with a parameter identification technique. For transient opera- 
tion, with the rotor controlled by the gyro, considerably larger hub moments 
were calculated than measured while calculated gyro angles were in reasonable 
agreement with measured gyro angles. Application of the parameter identifica- 
tion technique resulted in improved hub moment matching when as many as 10 
parameters of the model were adjusted simultaneously. The sensitivity of the 
blade response to small changes in the parameters was also calculated. The 
most sensitive parameters were found to be collective pitch, blade twist, tip 
loss factor, blade to gyro coupling constant. Lock number, and equivalent 
flapping hinge spring constant. 


INTRODUCTION 


The application of parameter identification techniques to the problem of 
estimating stability derivatives for fixed-wing aircraft has been reported by 
many investigators (refs. 1, 2, and 3), These techniques have been applied 
to helicopters more recently (refs. U and 5), In reference k 9 stability and 
control derivatives were estimated from helicopter flight data; in reference 
5, the feasibility of estimating specific parameters contained in a rotor 
math model was investigated. The results of reference 5 indicate that good 
estimates of the parameters can be obtained with identification techniques. 
However, in reference 5, the math model itself was. used to generate the data 
used in the identification procedure. This provides a check on the parameter 
identification algorithm, but provides no information concerning the validity 



of the math model itself. This limitation is removed in this report by use 
of data obtained in the wind tunnel for the identification procedure. 

A gyro-controlled hingeless helicopter rotor has been tested in the Ames 
4(3- by 80-Foot Wind Tunnel. The rotor blades were relatively stiff with an 
undamped natural flapping frequency of about 1.25 cycles /rev at 100 percent 
RPM. An expository analysis of this rotor system is given in references 6 
and T where calculated cyclic angles for trimming rotor hub and swashplate 
moments were compared with experimental data at high advance ratios (> 0.5). 
Vibratory loads on the same rotor at high advance ratios were studied in ref- 
erence 8. The mathematical model used iq the above three references contains 
simplifying assumptions that are valid only for high advance ratios. The 
mathematical model developed here is valid for advance ratios below 0.5. 
Although the wind-tunnel tests were aimed at the slowed-stopped operation of 
the rotor and therefore emphasized high advance ratios, some data were 
obtained for advance ratios less than 0.5. These latter data were used for 
parameter identification in refining the parameters of the math model of the 
rotor. In the math model, each blade of the rotor was allowed flapping and 
feathering degrees of freedom independent of the others. The first blade 
flapping mode was represented by rigid blade motion about a spring-restraint 
virtual hinge. The blades were also swept forward by an angle from the 

feathering axis, allowing flapping moments to be fed back into the feathering 
axis and hence to the control gyro. 

The behavior of the mathematical model for different inputs was investi- 
gated by integrating the time-varying periodic differential equations forward 
for several rotor revolutions. The simulation was compared with wind-tunnel 
data during both steady-state and transient conditions. Steady-state and 
transient responses for advance ratios up to 0.^0 were used in the comparison 
Parameters of the model such as blade inertia, flapping hinge spring constant 
blade to gyro coupling constant. Lock number, tip loss factor, collective 
pitch, and blade sweep angle were systematically varied using a parameter 
identification technique to give the best match between the simulation and 
the wind-tunnel measurements. The sensitivity of these parameters on the 
output response was also calculated. 


DEVELOPMENT OF THE ANALYTICAL MODEL 


An analytical model was developed for a three-bladed hingeless rotor 
with a high-speed gyro in the feathering system (fig. l). A list of assump- 
tions and the derivation of the equations are given in appendices A and B. 
The complete rotor equations are derived as follows. First, the equations 
for the blade motion are obtained; the equations for the gyro and gyro-blade 
interaction are then derived; and, finally, these two sets of equations are 
combined to yield the complete set. 


2 



Zs 




FWD 




Figure 1.- Gyro- controlled 
hingeless rotor. 


Figure 2,- Rotor-blade coordinate 
system. 


Blade Equations 


The equations for the flapping and feathering motion are obtained by- 
writing the moment equilibrium about the hub. Figure 2 shows the coordinate 
system used in writing the moment equilibrium. The X S ,Y S ,Z S coordinates 
represent the stationary axis system with the X s axis rearward over the 
tail section of the helicopter, the Y s axis to the right, and the Z s axis 
upward. The positive direction of rotation about each stationary axis is 
indicated in the figure. The azimuth angle rp is measured positive from the 
stationary axis X s in the direction of rotor rotation. The positive direc- 
tion of the flapping angle. 8, inplane angle C, and feathering angle 6 is 
also shown in figure 2. The effects of the inplane motion are included in 
the equations derived in appendix A. However, for the verification and dis- 
cussions in this report, the equations are simplified by fixing the inplane 
deflection at C Q , the preset blade sweep angle. This simplification is 
reasonable if resonance of the £ motion is avoided (ref. 6), 'The moments 
acting about the hinge point are the inertial, aerodynamic, and mechanical 
moments. Each of these moments for the flapping degree of freedom is devel- 
oped in this section. For the feathering degree of freedom, the equations 
are developed in conjunction with the gyro equations in appendices A and B 
because the feathering axis is mechanically connected to the gyro system. 

The inertia moment about the flapping axis is derived in appendix A 
(eq. (A13)): 




inertia 


= 18 

y 


- (i 


+ i ^ 
y o 


I }ft 2 8 
z 






-(i + i - i )qb - (i - i )c ^ 2 e (l) 

x y z z y o 
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The aerodynamic moment about the flapping axis is also derived in appendix A. 
For linear aerodynamics with no stall effects included, the aerodynamic 
moment is written in abbreviated form (eq. (A36)) as 


*V, 


aero 


as -Kl{-C2(^)3 + [-C1(U0 + C4(ip)tan 6 3 ] 3 + C3(ip)9i + C4{ip)0 

o 


+c4(t)0 + C5U)X + C6(ip)x + CT(^)X > 

eye o x y 


( 2 ) 


The coefficients of equation (2), defined by equations (A34) and (A35)> 
repeated here for completeness: 


K1 = 


yl tt 2 

L z. 

2 


are 


32 = ^1 


B3 

B4 

B5 


3 

B 4 

T 

5 



Cl(ifi) = £ q B 4 + p(cos ip + C Q sin ip )B3 + y 2 (c o + ? 0 cos 2ip + sin 2ip)B2 

C2(t|>) “ ^ + q (C Q Cos \p + sin tp)B3 

C3(i|>) = B5 + 2y(c o cos ip + sin ip)B4 + U 2 ^* - ^ cos 2ip + C Q Sin 2i^B3 

C4(*) = B4 + 2vi(£ q cos ip + sin ip)B3 + ^ cos 2ip + C Q sin 2 tAb2 


C5(ip) = B3 + y(C Q cos ip + sin ip)B2 


C6(ip) = b 4 cos ip + y(cos ip sin ip + C cos 2 ip)B3 

o 

C7(^) = b 4 sin ip + y(sin 2 ip + c, cos ip sin ip)B3 

o 


The reverse flow region encountered hy the retreating blade during for- 
ward flight was modeled in the aerodynamic moment term of equation (2). The 
reverse flow region was approximated hy assuming that reverse flow existed 
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from an azimuth 225° to 315° and from the root x * y. The appropriate 
changes required in the coefficients of the aerodynamic moment equation (2) 
are given in appendix A* 

The other moment component is produced by the spring restraint of the 
blade at the hub. For the flapping axis, this moment is given by k^(3-3 p ), 

where k represents the spring constant of the blade and (B-B ) represents 
p P 

the flapping deflection away from the precone position. 

The mechanical spring moment for the flapping axis can be written as 
(eq. (A37)): 


Vech = ' 6 p ) C 

The moment equilibrium about the flapping axis is now given by the sum 
of the inertia, aerodynamic, and mechanical moments. The flapping equation 
may be written by combining equations (l), (2), and (3) to form: 


^Vinertia + ^Vaero + ^"mech ^ 


which gives 


I 5 - (I + I C 2 ~ I )fl 2 B + | I 

y X y o z 


0 * D — ft 2 3 - (i + i - i )ae - (i - i k ft 2 0 

2 y F^-e x y z z y o 


-KL{-C2 (i/j ) 3 + [-Cl(*) + C4(i/»)tan <5 3 ]B + 03(1*00! + CU(xp)0 + Cb (\ *>)e 

O C- 

+C5(i|»)X + C6(ijOX + C?(ifOX } + k ( 3 - 3 ) = 0 
o x y 3 p 

(U) 


Equation (4) describes the flapping motion of each blade and will be used 
later to develop equations for a three-bladed rotor. 


Gyro Equations 

The gyro of the rotor system, modeled prescribes a feathering sequence to 
the blades. The equations for the gyro are developed about its axes system 
and the relationship of the feathering axes with the gyro is then established. 
The coordinates used to derive the gyro equations are shown in figure 3. The 
coordinates X S ,Y S ,Z S and X r ,Y r ,Z r are the stationary and rotating coordi- 
nates of the blade as previously defined. The X C ,Y C ,Z C coordinates define 
the rotating coordinates of the gyro. The equations for the gyro are 
obtained by expanding the angular moment tun equation [dH^/dt]^ = M Q . This 

equation is expanded in appendix B. In terms of the gyro coordinates 6j 
and S 2 , the equations for the gyro are (eqs. (B25) and (B26)): 
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I G «1 + (Xs" 2 ^ - I g « 2 )«1 + - 2ai G )6 2 - fiC g « 2 + C^j + _ M exCxc 

= V 1 sin(l|j+^) - V2COS (^+0 Q ) 

(5) 

I g «2 + (Ia^kp - X q S! 2 )«2 - ( 13 ^ - 2!2I g )«! + !JC s «i + C s S 2 + K g « 2 - M exCyo 

= vicos(i|;+^ o ) + V2sin(^+^ o ) 

( 6 ) 


The left sides of the equations 
represent the inertial properties of 
the gyro, the spring and damper effects 
(K s ,C b ) about the 6 ; and 62 axes, and 
the excitation of the gyro by the^ 
blades (M exCxc and M exCyc ). The right 

side of the equations represent the 
roll (v^) and pitch (y^ ) moment inputs 
of the pilot. The pilot input is via 
a moment (not displacement) applied to 
the gyro in the nonrotating frame. 

The gyro responds to the input moment 
as well as to moments applied by the 
blades. The M exCxc and M exc _^ c terms 

in equations ( 5 ) and ( 6 ) are the 
moments applied by the blades and must 
now be obtained from the feathering 
equations of the blade expressed in the 
gyro axes system. The moments contributing to the excitation are the blade 
inertia and mechanical moments about the feathering axis. Aerodynamic moments 
about the feathering axis are assumed negligible. In addition, flapping 
moments coupling into the feathering axis due to the forward sweep of the 
blade must also be included. The complete moment expression in the blade 
frame, given by equation (B27 ) , is expanded by substituting the appropriate 
moments and then expressing the moment equation in the gyro axis system. The 
excitation of the gyro by the blade expanded for a three-bladed rotor and 
transformed to the X C ,Y C coordinates is given as (eqs. (B3l) and (B32)): 

Mexc xo = C 2 {-I X ( § 1- I I - o e (&l- \ h- f S3) - (I z -x y )2 2 (a 1 - | e 2 

- | e 3 ) - k 6 ( ei - I 02- | e,) - (vyi.Jn^- § b 2 - | s 3 ) + d a -i y ) 

xc o fi 2 (gi- | 62- | e 3 )+C o k e [(Bi-Bp)- |(B 2 -B p )- |d3-B p )]J (T) 



Figure 3.- Gyro coordinates. 
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Mexo yo = ^ c 2 {-i x (e 2 -& 3 ) - c e ( 0 2 -e 3 ) - (i z -i y )n 2 ( 0 2 -e 3 ) - k 0 (e 2 -a 3 ) 

-d +i -i )n(e 2 -6 3 ) + (i -i )? n 2 (6 2 -s 3 ) + c k [(b 2 -b ) - (e 3 -e )]} 

x y z zyo op p p 

( 8 ) 


The gyro equations are now obtained by combining .equations ( 7 ) and (8) with 
(5) and (6) and transforming the feathering angle 0 for each blade and their 
derivatives to gyro coordinates. The transformation between 0 j, 0 2 j an< 3 - 
6l, 62 is given by (eq. (A 4 o)): 

©i = C2 icos(i-l) — g* + 5 2 sin(i-l) ; i = 1 , . . b ( 9 ) 

where b = 3 . Substituting the above transformation and its derivatives into 
equations (T) and (8) gives the complete gyro equations for a three-bladed . 
rotor (eqs. (B 35 ) and (B36)): 

(y- I c|ljs i+ | 0 e cl6 1 +C s 5 1 +(X 3 n 2 k R -I G n 2 )6 1 + | c|(l z -I y )B 2 6 1+ fkgCf^+K^i 

+ (l 3 flk B -2!lI G ,)« 2 -n0 8 6 2 +C 2 (l x +I y -I z )n^i- | B 2 - § B 3 )-C 2 (l z -I y )c 0 !2 2 

x(B t - I 3 a- | 6a)-C2C 0 k e [(Bi-B p )- § ( 3 2 - 3 p )- | ( 33 - 3 p )] 

= Visin(ifj+ip o )-v 2 coB(i(J+^ o ) 

( 10 ) 

( I G + I C i I x>' 2+ I <= 8 cli2 + C s «2+(l 3 !2 2 k E -R 2 I G )S2+ | cl(l a -I y )n 2 6 2 + |k 6 ci« 2 +K s 5 2 

-(l30k R -2flI G )6 1 +S2C g {i+ ^ C 2 (l x +I -I z )n(@ 2 -0 3 )- ^ C 2 (l z -I y )C o !J 2 (S 2 -B 3 ) 

y/o" 

- C 2 c k [(e 2 -$ ) — ( 3 3— 3 )] = v 1 cos(^+Tp )+v 2 sin(\p+ip ) . 
d o p p p o o 

( 11 ) 


Rotor Equations 

The flapping equation ( 4 ) was written for one blade. For the three- 
bladed rotor, the flapping angle 3 and the feathering angle are subscripted 
1, 2, and 3; for the azimuth angle tp in equation ( 4 ), the substitution of 
<p s ip+ 120 ° , and ip+ 240 ° for each blade must be made, respectively. The flapping 
equation can then be written as 
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i g -(i +i 5 2 -x )sj 2 g. + 1 1 5 s - n 2 e.-(i +i -x ) 08 .-(i -i ); fi 2 e. 

y i x y o z i 2 y R-e 1 xyz 1 zyo i 

-KlJ-C2[i|j+(i-l) ^Js-Cl[ip+(i-l) ^.]3 i +C4[^+(i-l) ^-j^tan 6 3 +C3[V(i-l) 

+ cu[^(i-i) f] (e o + V )+C5 ^ (i_1) f] X 0 +C6 t +(i -l) f] X x +CT [* +(i “ l) f] X y) 

+kJg.-B ) = 0 i = 1, 2, 3 (12) 

PIP 


The feathering angle 6 and its derivatives in equation (12) are written in 
gyro coordinates (see eq. (9)). In appendix A, the transformation equation 
(9) has been applied to equation (12) to convert the feathering angles to gyro 
angles. The equations for each blade are given in expanded form by equations 
( A4l ) , (a 42), and (a 43). The gyro equations with the feathering equations 
interconnected for the three blades were written in gyro coordinates in the 
previous section. The complete set of rotor equations is now given by the 
two gyro equations (10) and (ll) and the blade equations from equation (12) 
for i = 1, 2, and 3. 


Method of Solution and Transformation of Equations 
to Measurement Coordinates 


In this section, the transformation to measurement coordinates is dis- 
cussed, followed by a description of the method for calculating time histories 
from the math model. 


The five equations discussed in the previous section (eqs. (10) and (ll) 
and eq. (12) for i = 1, 2, 3) are next written in first-order form for the 
digital simulation. The first-order state variables are defined as follows: 


— -1 
*1 


■ 

6 i 

X 2 


Si 

x 3 


62 

x 4 


6 2 

x 5 


Pi 

X 6 


61 

X? 


32 

X 8 


32 

Xg 


33 

x 10 


h 

L J 


L J 


(13) 
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Applying these transformations of state variables to the five rotor equations 
gives the equations of the rotor in the form x = A(t)x + ?(t) (see table l). 
The equations were integrated using the Runge-Kutta, fixed stepwise integra- 
tion routine. The output of the integration process produced flapping angles 
and gyro angles. However, flapping angles could not be measured during the 
wind-tunnel tests; hub moments were recorded instead. Therefore, a transfor- 
mation was required to convert the flapping angles to hub moments. The trans- 
formation from moments to flapping angles, derived in appendix C (eqs. (C12) 
and (C13)) 9 is given below: 

m b x = k e(i ^ e 2 - | *'3. S 3 ) d 1 *) 



where the shearing force at the hinge was considered small since the rotor 
vas lightly loaded. Equations (lU) and (15) were written in the form 
y = Sx, where y is the vector of moments and^ x is the vector of flapping 
angles. The complete set of rotor equations (x = K( t)x + f(t) and y - Cx) in 
first-order form is shown in table 1. 

The simulation and measurement results were both performed in the rotat- 
ing coordinate frame X r ,Y r ,Z r (shown in fig. 2). During the analysis of the 
data, it was found in some instances that the hub moments and angles could be 
better understood if the moments and angles were transformed to a nonrotating 
coordinate frame at the hub center - the X S ,Y S ,Z S frame discussed previously 
(fig. 2). The transformations between these reference frames are also given 
in appendix C. 


BEHAVIOR OF SIMULATION MODEL 


Parameters of Model 


The nominal values of the parameters for the rotor and the gyro used in 
the analytical model simulation are given in table 2. These parameters were 
obtained by direct measurement or calculation. 


The first flapping mode of the blade was represented by rigid blade 
motion about a spring restraint virtual hinge located a distance e radially ■ 
from the hub center. The spring constant about this virtual hinge was calcu- 
lated from 

k e = V 2 ( p2 - 1 - 1 rr? - 8 tan 4 3) (16) 

where ' 


I 


y 


blade flapping inertia 


2 rotor angular velocity 
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p flapping natural frequency to rotor frequency 

e virtual hinge location 

R. "blade radius 

y Lock number 

63 pitch flap coupling angle 

However, for this rotor, 63 was zero and the spring constant was evaluated 
from the first three terms of equation. (l6). The values for p in equation 
(l6) were obtained from reference 6 and are given as a function of percent 
rotor RPM in figure h. The values for I y , fl, and were obtained by mea- 
surement, while the virtual hinge location e was found to be between 11 to 
15 percent of the blade radius (to be discussed later). The variation of 
spring constant with percent rotor RPM with Iy = 363.^ kg-m 2 (268 slug-ft 2 ), 
63 = 0, and e/R^ = 0.11 is also shown in figure U. 



percent, rpm 


Figure U.- Blade flapping frequency 
ratio and blade spring constant 
as a function of percent rotor 
RPM (100-percent RPM = 355.6). 


Time Response 

The rotor flapping response for 
different inputs was investigated by 
integrating the time-varying periodic 
differential equations forward for 
several rotor revolutions. The param- 
eters in the math model used to gen- 
erate these time histories are given 
in table 2 and those unique to the par- 
ticular case studied are given with 
the results in each figure. The vir- 
tual hinge location was fixed at 11 
percent of the blade radius to illus- 
trate the model behavior. First, the 
response of the rotor in hover was 
obtained. Figure 5(a) shows the blade 
flapping angles before and after a 
step input of 2° was applied to the 
collective control at the start of the 
second rotor revolution. The gyro 
angles and the hub moments are zero 
since the swashplate moment inputs 


were taken as zero for this run. The transient motion of the blade flapping 
angles is shown for about one rotor revolution (fig. 5(a)) after which the 
coning angle reaches a new steady-state value. Figure 5(b) shows the 


response of the rotor in hover but with an incremental step swashplate roll 
moment input of = 2T Nm (20 ft-lb) applied at the second rotor revolu- 

tion (Ma = hi Um (30 ft-lb) before the step input). The swashplate pitch 
■ sp 

moment and the collective pitch angle were held fixed at Ma = -l4 Nm 


10 



(-10 ft-lb) and 0 Q = 8.57° * respectively. For this case, with moments 
applied to the gyro, the hub moments are not zero. The transient response of 
hub moments in rotating coordinates (R- frame), gyro angles, and blade flapping 
angles are shown in figure 5(b). 

Two additional transient responses were calculated but for the rotor in 

forward flight with advance ratios of u = 0.212 and 0.485 (figs. 6 and 7). 

Figure 6 shows the transient response of the blade flapping angles, gyro 

angles, and hub moments for a step moment input applied about the swashplate 

roll axis. The first two rotor revolutions in figure 6 show the steady-state 

behavior with = 4l Nm (30 ft-lb), Ma = -lU Nm (-10 ft-lb), and 

Y sp Y sp 

0 Q s 8,57°. At the second rotor revolution, an incremental step input about 
the swashplate roll axis of 27 Urn (20 ft-lb) was applied while the swashplate 
pitch moment and collective pitch angle remained the same. The transient 
motion that followed is shown from the second to the eighth rotor revolution. 




(a) Response to a step collective pitch input at the start of the 

second rotor revolution. 

Figure 5.- Transient response at hover. 
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(b) Response to a swashplate roll moment input at the start 
of the third rotor revolution. 

Figure 5.- Concluded. 


In like manner, the same swashplate inputs were applied as when the rotor was 
operating at an advance ratio of y = 0.1*85 (fig. 7). A comparison of the 
transient behavior of the rotor at the two different advance ratios indicates 
that for the lower advance ratio the response is sinusoidal, whereas for the 
higher advance ratio the response shows the presence of higher harmonic terms. 
This effect is due to the periodic terms in the differential equations which 
become significant at the higher advance ratios. 
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Rotating coordinates - hub moments 
Pitch, ft -lb Roll, ft-Jb 




^.Z\Z 

k£ = 184394 Nm/rad 
(136p02 ft-lb/rad) 
y- 4.74 

V = 36 m /sec (70 knots} 

£1 = 33.75 rad/sec 
(90.6% rpm) 

C fl = 9.5Nm/rod/sec 
(7ft-lb/rad/sec) 

C, -)08.5Nm/rod/sec 
(aOft-lb/rad/sec) 
fcU = 40.7 Nm to 67. 8 Nm 
p (30 fl -lb to 50ft-lbl 
* -l3.6Nm(-IOft~lb) 

s p 

e/R b = .H 



Figure 6.- Transient response to a swashplate roll moment input at the start 
of the third rotor revolution for u = 0.212. 


13 


6000 



» - 6000 L 




-6000 


J. 


J 1 1 I I I 



ft - .485 

kg - J38p6l Nm/rad 
(101.^43 f 1-lb /rad) 
y = 4.57 

V = 46,3 m/sec (90 knots) 
12 ~ 18.96 rod /sec 
(51% rpm) 

C g = 9 . 5 Nm /rad/sec 
(7ft-lb/rad/5ec) 

C s =108.5 Nm/rod/sec 
(80ff-lb/rod/sec) 

M,* = 40.7 Nm 10 67. 8 Nm 

p (30 ft -lb to 50fMb) 
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Figure 7 .- Transient response to a swashplate roll moment input at the start 
of the third rotor revolution for y = O.H85. 


APPLICATION OF PARAMETER IDENTIFICATION TECHNIQUE 


The parameter identification process is useful in redesigning and refin- 
ing the parameters of the math model using wind-tunnel or flight measurements . 
The parameters to he determined are those that appear in the coefficients of 
the differential equations, such as inertia, spring constant. Lock number, 
and damping constant. The general outline of the identification process is 
sketched in figure 8. The same input u is applied to "both the math model 
and the plant ("plant” refers to the helicopter rotor test in the wind tun- 
nel). The best values of the parameters of the math model are those that 
minimize the difference between the model and plant outputs {f and y). The 
identification process is an algorithm that compares the outputs of the math 
model and plant and adjusts the parameter values (inertia, spring constant, 

Ik 


damping, etc.) until this difference 
is a minimum. What this minimum 
should be is established by least- 
squares fit within the limits of 
table 2, and the identification algo- 
rithm iterates until this criterion 
is reached. It may not always be pos- 
sible to reach this minimum or even 
to obtain an improvement in the match- 
ing of the model and plant outputs. 
This report does not address these 
problems, but will make use of the 
technique in obtaining, where possi- 
ble, a better match between the simu- 
lation and the measured data. For a 
detailed discussion of the parameter 
identification theory as applied here, 
see reference 9. 

The identification program had the capability of identifying as many as 
36 parameters of the mathematical model. For the results that follow, 6, 7, 
and 10 parameter cases were studied. The choice of parameters was based on 
the accuracy to which the parameters were known, the sensitivity of the param- 
eters in affecting the response, and whether the parameters were linearly 
independent. The accuracy to which the parameters were known a priori estab- 
lished lower and upper bounds of variation for each parameter. A sensitivity 
analysis of all parameters indicated which ones were most sensitive in affect- 
ing the response. The linear independence or uniqueness of the parameters cal- 
culated by the identification program is discussed in detail in reference 9. 
Thus the parameters studied here were chosen on the basis of a priori knowl- 
edge, sensitivity, and linear independence. Tables 3 and k present the param- 
eters identified and their nominal values obtained from reference 6. 

The computation time required to identify more than 10 parameters over 
several rotor revolutions is large. The approximate time required to estimate 
one parameter using one measured data point per iteration was found to be 
about 1/20 of a second on an IBM 360/70 computer. This time could vary con- 
siderably since the computation time is not a linear relationship of the 
parameters estimated. Adequate reconstruction of the measured waveforms for 
parameter identification purposes could be achieved by sampling the data every 
15° in azimuth. The integration step size was allowed to become no larger 
than 2°. In some cases, a larger step size could be used, but the 2° step size 
was found to be satisfactory in eliminating significant error accumulation 
when the equations are integrated for several rotor revolutions. 


New values 



Figure 8.- Block diagram of parameter 
identification technique. 
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COMPARISON OF MEASURED AND SIMULATED DATA 


Discussion of Measured Data 

The measured data used to check the analytical model were obtained in 
the Ames U0- "by 80-Foot Wind Tunnel for a three-bladed, hingeless rotor with 
a gyro coupled to the blade feathering axis. Figure 9 is a photograph of the 
model in the wind-tunnel test section. This model was used to investigate 
the behavior of the rotor under slowed-stopped conditions. For these condi- 
tions, the rotor was operated near zero lift. Hence the measurements taken 
during tunnel testing are for a lightly loaded rotor at relatively high 
advance ratios. The following measurements were recorded: load cell moments 

(moments in stationary coordinates 0.91^ m below the hub), shaft bending 
moments , swashplate moments and angles , and blade angles . Data were available 
to check the analytical model and to identify parameters at forward speeds 
from about 31 to Ul m/sec (6l to 80 knots) and rotor RPM from l48 to 197. A 
complete description of the test and data collected is given in reference 6. 
For the steady-state cases used to check the analytical model, the rotor was 
operated in a locked gyro mode. For this mode, the swashplate was held fixed 
at specified angles while the rotor flapping angles were allowed to come to a 
steady-state condition. The transient data had been recorded by taking the 



Figure 9 . - Experimental rotor rig in the Ames U0— by 80-Foot Wind Tunnel. 
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rotor out of trim and then allowing it to return to a trim condition. This 
procedure amounted to a step change in the swashplate angle. The purpose of 
this procedure was to establish the rotor transient independent of gyro 
motion. 

After the measured data were digitized and converted to engineering 
units, it became apparent that some filtering of the data was necessary, A 
low-pass digital filter with a corner frequency of 1 6 Hz was used to process 
the transient data. The low-pass filter eliminated some of the high-frequency 
noise and vibrations from the measured data without introducing any noticeable 
phase shift for frequencies of interest here. Figure 10 compares' the filtered 
and un filtered data for a typical transient run to show the quality of the 
data and the results of filtering. 

In comparing the digital simulation with the measured data, care had to 
be taken to use the correct initial conditions in the model simulation. The 
initial conditions required for the computer simulation were gyro and flap- 
ping angles and rates and the moments applied to the swashplate. Since only 
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Figure 10.- Typical measured transient data with and without filtering; 

y = 0.355* V 36 m/sec (70 knots). 
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gyro angles and swashplate moments had been recorded for each record taken in 
the -wind tunnel* the gyro rates and. flapping angles and rates had to he 
established. The gyro rates and the flapping angles and rates' were obtained 
by entering into the computer simulation the known values of swashplate 
angles, rates, and moments and integrating the equations to steady-state con- 
ditions. The flapping angles and rates thus obtained served as the initial 
conditions for the particular .tunnel record compared. 

The next section presents steady-state and .transient runs for different 
forward velocities and results obtained by applying the parameter identifica- 
tion program to the measured wind-tunnel data. For a steady-state condition, 
the gyro was held at a fixed attitude (locked gyro mode). For the transient 
run, the gyro was free to move but under the influence of moment inputs to 
the swashplate. 


Steady-State Case: Locked Gyro 

Figures 11 and 12 compare measured (discrete points) and simulated 
(solid line) hub moments for two different forward velocities but with about 
the same advance ratio. The main difference between these two cases is in p 
the dimensionless flapping frequency. The rotor was operated in a locked 
gyro mode while the measured hub moments were recorded. In like manner, the 
mathematical model simulation was forced with the gyro angles and moments 
recorded during tunnel testing while the hub moments were allowed to vary. 

The nominal values used in the simulation are given in figures 11 and 12 for 
a virtual hinge e located at 11 percent of the blade radius. The math 
model shows similar waveforms for the hub moments without the presence of 
higher harmonics in the measured data. The vehicle moments are the 1/rev 
moments in rotating coordinates and these moments are important in control- 
ling the vehicle. The higher harmonics seen in the measured data are asso- 
ciated with' the vibratory modes of the vehicle. An amplitude difference in 
the math model results is apparent for both advance ratios studied. The 
largest difference in amplitude occurs for lower forward velocity with 
U = 0.402 and p = 1.62 (fig. 12). 

The parameter identification program was applied and a number of param- 
eters were adjusted simultaneously to improve the hub moment match. The 
results of the hub moment matching obtained after 10 parameters of the math 
model were identified are shown in figures 13 and l4, with the nominal values 
and estimated parameters for each case shown in tabular form. The inflow, 
uniform for both cases, was calculated from equation (A3l). The spring con- 
stant of the blade was decreased from the nominal value. The ratio of hinge 
offset to blade radius (e/R^) was allowed to vary and was found to be at 11 
percent of the blade radius. The remaining parameters varied from 2 to 3 
percent, which was the accuracy to which these parameters were known a priori 
Improvement in the huh moment match is shown in figure 13 after the param- 
eters were varied from their nominal value. The large difference between the 
measured and math model results for y = 0.402 (fig, 12) was considerably 
reduced as shown in figure l4. A comparison of the identified parameters, 
for the two steady-state cases studied, shows similar values of those 
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Rotating coordinates-hub moment 


Nominal values of parometers 
V = 41.2 m/sec (80 knots) 


rpm = 197 

1 20.66 rod/sec (55.5% r^m) 
f> i.l83kg/m3(.00229^l-) 

f f 0 

C T = .00097 

- 141,092-^(104,064^—^) 
rod rad 

y = 4 67 
*S P = 2.15“ 

8o = -3. 50° 

5 p 

= 98.3 Nm (72.5 ft-lb) 

Sp 

M 8 - -34 2 Nm (-25.25 ft-lb) 

Sp 

e/R b * .11 




Nominal values of parameters 
V- 31.3 m/sec ( 61 knots) 
rpm * 148 . 

ii = 15.50 rad/sec (41.6% rpm) 
p - 1.188 kg/m3 (.00231 slugs/ftS) 
C T = .0015 

1*0= 125,214 Nm/rad (92?53 ft-lb/rad) 
y = 4.69 
<£sp= 3.42° 

8 S p=-66 0 

M^, sp = 45.8 Nm ( 33.8 ft-lb) 

M«sp' 57 4 Nm (42.3 ft-lb) 
e/R b = .11 



Measured data 

— Moth model 
using nominal 
values from 
table 3 



Figure 11.- Measured and simulated Figure 12.- Measured and simu- 

hiib moments in rotating coordi- lated hut moments in rotating 

nates; u = 0.399, p = 1.45. coordinates; v - 0.402, 

p = 1 . 62 . 
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Nominal value of parameters 


identified parameters 
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8.57° 
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363.9 kg-m 2 (268.4 slug- ft 2 ) 


4.67 


ff 0 = 7.97° 
B j, -- .95 


C 2 = .85 
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l y = 360.2kg- rT! 2 (265. 7 slug- ft 2 ) 

Y = 5.13 


I4IP92^ (104,064^1 
1.5° 

2.25° 

.tl 


V 129,345 ^(95f 00^) 
£o 8 1-6“ 

2.36” 
e/R b = .11 



— — — Measured data 

Math model 

using identified 
parameters 


Figure 13.- Measured and identified hub moments - 
uniform inflow; u - 0.399, p = 1*^5. 

parameters allowed to vary. The inertia, tip loss factor. Lock number, blade 
sweep angle, collective pitch, gyro coupling constant, and hinge offset were 
the same after the identification for the two steady-state cases studied. 


Transient Case: Free Gyro 

In this section, the math model simulation output is compared with a 
transient run recorded during wind-tunnel testing. In addition, the param- 
eter identification algorithm was again applied to improve the matching of 
simulated and measured data. 

The transient measurements were recorded with the rotor operating in a 
free gyro mode and with a step roll moment applied about the swashplate roll 
axis. For the free gyro mode of operation, the gyro is free to move in 
response to moments from the blades and swashplate. The step input applied 
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to the roll axis of the swashplate 
was accomplished by taking the rotor 
out of trim and then moving the 
swashplate roll moment control 
quickly back to a trim position. 

The transient run was obtained at a 
forward velocity of 36 m/sec (70 
knots) with a rotor speed of 192 RPM. 
The measured data had been filtered 
(discussed earlier) to eliminate the 
high-frequency components found on 
some of the recorded data. The nom- 
inal values of the parameters used 
in the math model simulation are 
shown in table 2. Those parameters 
unique for this advance ratio are 
tabulated in the figures next to the 
measured response. 

The measured response (discrete 
points) of swashplate roll and pitch 
moments , hub moments , and gyro 
angles is shown in figure 15. The 
mathematical model response (solid 
line) of hub moments and gyro angles 
is superimposed on the measured data. 
Both measured and calculated hub 
moments and gyro angles are presented 
in rotating coordinates. The swash- 
plate roll and pitch moment inputs 
used as forcing functions in the 
mathematical model simulation are 
shown in the upper two traces of 
figure 15. The step change of the swashplate roll moment is seen at the end 
of the third rotor revolution while the swashplate pitch moment remained 
unchanged. The initial conditions were obtained by forcing the math model 
with the swashplate inputs prior to the step change in the swashplate roll 
moment. Integrating the equations with these swashplate inputs until steady 
state was reached produced the initial conditions required to start the math 
model simulation. With these initial conditions, the measured and simulated 
time responses for this step swashplate roll moment input are shown in the 
lower four records of figure 15. Note that the gyro angles from the math 
model simulation match the measured data fairly closely. However, the magni- 
tude of the hub moments is considerably different for the simulated data com- 
pared to the measured data. 

Several identification runs were made so that a different set of param- 
eters could be identified to improve moment and angle matching. Figure l6 
shows the matching of hub moments and angles when six parameters of the math 
model were allowed to vary. Each parameter was allowed to vary over a per- 
centage range from their nominal values. The percentage variation was 


Nominal values of parameters 
8 0 -8.57° 

B<=.97 
Cg *.87 

I j, *363 4 kg-m 2 (268s!ug-ft z ) 
I z * 363.9 kg- mZ< 268.4 slug -ft 2 ) 
y - 4.69 

K * ='23214 7§T < 92 ? 53 75? » 


C„ =1.5° 

% = 2.2 5° 
e/Rb = . 1 1 



Identified parameters 
8 0 *7.97° 

Bf = .95 
C g = .85 

1 y = 359.7 kg-m z (265.3 slug -ft 2 ) 
I z = 360.2 kg-m z (265.7slug-ft 2 ) 
y = 5.16 

k fj - 8 7 651 (64646 ) 

Co =1.6° 

$ p *2.14° 
e/Rb * . II 


Measured data 

Math model 

using nominal 
values from 
table 3 


2000 


-2000 



Figure lU.- Measured and identi- 
fied hub moments - uniform 
inflow; u = 0.U02, p = 1 . 62 . 
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Nominal values of parameters 
I* =.54kg-m2(.40slug-ft2) 
I y = 363.4 kg-m 2 
(268 slug-ft 2 ) 

I z =363.9 kg -m2 
(268.4 slug-ft 2 ) 

I 3* .41 kg-m2(.3slug-fl2) 
Iq= ,2lkg-m2(.l5slug-fl2) 
k^ 138752 Nm/rad 
(102?38 ft-lb/raa) 
y- 4 57 

C s = 108.5 Nm/rad/sec 
( 80 ft-lb/rad/sec) 

Cg = 9.5 Nm/rad/sec 
{ 7 fMb/rcd/sec) 

Rp * 5.03m (16.5ft) 
e/R b = .11 

C= ,357m (1.17 f f ) 

B, = .97 

£0 = 15 ° 

ft = 20.l7rad/sec 
(542%rpm> 

83= 0° 

S| = -9.43° 

0p* 2.25° 
p = 1.158 kg/rr3 

(00225slug/ft 3 ) 

C 2 = 87 
9 0 = 0.57° 

Vo = 60 ° 

O rprT i= I0200rpm 
k e = 0 
k s = 0 

V* 36 m/sec (70knots) 



Figure 15.- Comparison of math model simulation with measured data, \i - 0.355. 
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Nominal parameters 
6 0 * 8.57° 

0, * -9.43° 

B/ * .97 
Y * 4.57 

- f38,752 Nm/rad (102^338 fHb/rad) 
So s 1-5° 


Estimated value of parameters 
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established from the accuracy with which these parameters were known a priori 
(table 2). Most parameters were known between 2 and 5 percent, although some 
(such as virtual hinge location and spring constant) were, known far less 
accurately and therefore were allowed to vary over a large range. The hub 
moment matching was considerably improved after only a few iterations, while 
the angle matching remained about the same. The phase relation of the hub 
moments between the calculated and measured data did not improve. The param- 
eters that gave a better match of the time histories are also shown in 

figure l6. 


Figure IT shows the matching obtained when seven parameters were allowed 
to vary. In this case, all parameters were the same as in the previous case 
except for the blade twist, which was held fixed while the flapping inertia 
and hinge offset were additional parameters that could be adjusted. The 
match of hub moments is not too different from the previous case. Slightly 
better phasing of calculated and measured data was obtained for the first 
three rotor revolutions* The parameters for this match are also shown in 

figure IT. 

One additional case was studied. This case differed from the previous 
one in that the simulation was forced with the same gyro angles as had been 
measured. The same parameters as for the seven-parameter case were allowed 
to vary; the precone angle, advance ratio, and the mechanical advantage con- 
stant were also allowed to vary from their nominal values. A total of 10 
parameters could be adjusted by the identification program. Figure 18 shows 
the comparison for this case and the corresponding identified parameters are 
given in the figure. . The gyro angles are in exact agreement in this case 
since the math model was forced with the measured gyro angles. The hub 
moment match is similar to that obtained in the previous case. The identifi- 
cation algorithm forced the magnitude of the hub moments to be the same, but 
better phasing of calculated and measured response over 10 rotor revolutions 
was not obtained. Limitations in the mathematical model or in the measured 
data used for the identification might have resulted in the lack of phase 
matching. 


Discussion of Results 

The parameters obtained from the identification program for the steady- 
state and transient cases studied are shown in tables 3 and U. The tabula- 
tion of the parameters for each advance ratio shows a consistency in the mag- 
nitude of the estimated parameters. For the two steady-state cases investi- 
gated, the estimated values of 10 parameter cases were similar. For the 
transient run studied, some variation in the magnitude of the parameters was 
found when a large number of parameters were estimated. Some of the scatter 
in the estimated values of the T and 10 parameter cases may have been due to 
noise present in the measured data. Also, for both the steady-state and 
transient runs, the rotor was not lifting, which could make the contamination 
of the measured signal by noise more significant* Two other factors that may 
have influenced the lack of matching are deficiencies in the mathematical 
model and the range over which the parameters were allowed to vary in the 
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Rotating coordinates -hub moments 
Pitch, ft-lb Roll, ft-lb 


Nominal parameters 

00 = 8.57 
B* =-97 

I y =363.4kg-m2(268slug-ft2) 
x = 4.57 


Estimated value of parameters 

0o*9.I7 
% = .985 

I y = 366.5 kg-m2 (270.3 slug-ft2 ) 
Y = 4.0 


k 0 = 138752 Nm/rad (I02?38ft-lb/rad k fl = 97,127 Nm/rad (71637 ft-lb/rad) 
Co= 1.5° So 8 1-8° 

e/Rtr .11 e/R b = .15 _ 


Math model simulation with 
estimated parameters 





Figure IT.- Comparison of math model simulation 

seven estimated parameters, p 


with measured data using 
= 0 . 355 . 
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Rotating coordinates - hub moments 
Pitch, ft — lb Roll, ft-lb 


Nominal parameters 


Estimated value of parameters 
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C 2 = .87 

I y =363.4 kg-m2{ 268 slufl-ft®) 
y = 4.57 


G 0 = 9.17 
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C 2 = .85 

I y =359.7 kg-m2 (265.3 slug-ff2) 
y = 5.02 


138,752 Nm/rad ( 102^38 ft-lb/rad) k 0 = 97,795 Nm/rad ( 72,130 ft-lb/rad) 


/j. = .355 

Co = 1.5° 
/3 p * 2.25° 
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ft ~ .348 
Co - 1-8° 
P P =2.I4° 
e/R b = .15 


Math model simulation with estimated 
parameters — fixed swashplate angles 
Measured data 


-4000 
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Figure l8.~ Comparison of math model simulation with measured data using 
10 parameters with fixed swashplate angles, y = 0.355. 



identification program. Only the flapping and feathering degrees of freedom 
were modeled here. But the inplane degree of freedom as well as blade tor- 
sional effects may be required in the mathematical model to obtain better 
matching of the measured data. The other factor influencing the time history 
matching is the limit imposed on the parameters over which the identification, 
program could search. The limits on the parameters were set by the accura- 
cies with which the parameters could be measured or estimated. An example is 
the spring constant of the blade, which is not constant for all rotor RPM and 
therefore is more difficult to determine. The virtual hinge location e is 
another parameter whose limits were not well defined. These factors will 
influence the time history matching obtained. 

Note again that the identification in all cases was carried out over the 
first 2 or 3 rotor revolutions. But the comparisons of the simulated and 
measured data using the identified parameters are presented for 10 rotor rev- 
olutions. This was done to show that an improved match for a very short por- 
tion of the measured data resulted in an improved match throughout the entire 
transient range of the measured data. Convergence of the parameters to be 
estimated was usually obtained after five iterations. Larger numbers of 
iterations were performed and larger sets of parameters were estimated, but 
hub moment and gyro angle matching were not improved. 

The identification program also calculates the sensitivity of the output 
to each parameter. The sensitivity was defined (ref. 9) as a dimensionless 
ratio of the change in the output due to a change in the parameter. For this 
model, the outputs are the roll and pitch hub moments and the roll and pitch 
gyro angles in rotating coordinates. A list of the sensitivities for the 
model parameters is given in table 5. Parameters of high sensitivity are 
collective angle (8 0 ), twist angle (6^ of the blade, tip loss factor (B^) , 
and blade to gyro coupling constant. The sensitivity of the next three param- 
eters in order of decreasing sensitivity are Lock number (y) (or blade 
inertia Iysl z )» spring constant (kj, and advance ratio (u). Accurate mea- 
surements of collective angle, twist angle, and blade flapping inertia can be 
obtained as well as the Lock number (which can be analytically evaluated). 

But the values of spring constant as modeled here are difficult to evaluate 
since they are not constant for different flight conditions. The best method 
for obtaining improved estimates for some of these parameters is through the 
use of parameter identification techniques. The information gained from the 
sensitivity analysis portion of the identification technique used here can be 
very useful in determining where improvements in measurement or the analyti- 
cal evaluation of parameters should be made. 

It was found from the application of this identification technique that, 
to obtain refinements of the parameters for the mathematical model, two con- 
ditions must be satisfied. First, the mathematical model of the physical 
system must be complete except for constant parameters. Second, the values 
and the range of uncertainties for each of the unknown parameters should be 
specified either from analytical techniques or from measurements. Subject to 
these conditions, the identification technique is applicable to the estima- 
tion of parameters of complex rotor systems from wind-tunnel or. flight data. 
Moreover, the utility of the method is evident in the identification of those 
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parameters that cannot he measured directly (such as the equivalent flapping 
hinge spring constant ) , hut are computed indirectly by. the algorithm through 
the complex relationship between the parameters and the measured response* 


CONCLUSIONS 


Comparisons of wind-tunnel data for advance ratios up to O.UO with the 
mathematical model simulation indicate that the mathematical model derived 
here gives a reasonably good representation of the rotor hub moments in 
steady-state operation. The match between the mathematical model and. the 
measured data was improved by adjusting selected model parameters with param- 
eter identification techniques. As many as 10 parameters were estimated 
simult aneously . 

For the transient runs with the rotor operating in a free gyro mode, con- 
siderably larger hub moments were calculated than measured. However,, gyro 
angles during the transient operation were in close agreement with the math 
model simulation. An improved match in the magnitude of the hub moments 
could be achieved using parameter identification techniques, but phase match- 
ing of the time histories over a large number of rotor revolutions could not 
be obtained. More extensive mathematical modeling of the rotor is necessary 
to further improve the time history matching. 

The sensitivity of the blade response to small changes in the parameters 
was calculated. The most sensitive parameters were collective pitch, blade 
twist, tip loss factor, blade to gyro coupling constant. Lock number, and 
equivalent flapping hinge spring constant. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9^035, June 5, 19TU 
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APPENDIX A 


DERIVATION OF BLADE EQUATION 


The rigid blade equations are derived in three parts: first, the 

moments due to inertial forces; second, the aerodynamic moments acting on th€ 
blade; and third, the mechanical moments that balance the inertial and aero- 
dynamic moments on the. blade. Figure 1 shows the gyro-controlled rotor to be 
modeled. 


Inertia Moments 

Figure 2 shows the right-handed coordinate system, used to derive the 
equations for the blade. The X S ,Y S ,Z S coordinates define the stationary 
coordinates of the blade. The X s axis points rearward over the tail sec- 
tion of the helicopter while the Y s axis points out of the right side of 
the hub. The Z s axis points upward following the completion of the right- 
handed coordinate system. The coordinates X r ,Y r ,Z r rotate at angular veloc 
ity Q of the rotor whose reference point Is taken from the stationary axis 
X s . The positive directions of the flapping angle 3, inplane c, and fea- 
thering angle 0 are defined in figure 2. 

In developing the inertia moments, the following definitions and assump- 
tions were made: 

(1) Only the first flapping mode is considered; it is represented by 
rigid blade motion about a spring-restrained virtual hinge located a distance 
e radially from the huh center. 

(2) Blade feathering (pitch) Is represented by rigid-body torsional 
motion. 

(3) All deflection angles (£„ 0, £) are small. 

(4) The inplane degree of freedom is not represented. 1 

The B frame (X^ ,Yt, 9 Z ^ ) is defined by the rotations of 3, and 0. 

The equations are developed in the B-frame since it represents the principal 
axes of the blade. The equations are then transferred to the R- frame (X r , 
Y r ,Z r ), which is also referred to as the shaft axis system of the blade. 

The rotational equations of motion in the B-frame (blade frame) are 
obtained by expanding the angular momentum equation for a rigid body (see 
ref. 10): 


M as I W + 

X XX 


(I - I )to to 
z y y z 


(Al 


M = I to 

y y y 




I )tl) U) 

z x z 


(A2 


where the effect of hinge offset e will he added later. To further expand 
these equations > the angular velocities and angular rates must he expressed 
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in terms of the coordinates shown in figure 2. Since only the feathering and 
flapping equations are developed here, the angnlar velocities ,co 2 and 

the angular accelerations and cy are needed. In developing these terms, 
small-angle approximations are made and terms small compared to 1 are dropped, 
which gives (ref. ll) 


w *.ne - - B + £e 

y 

co = n + 30 + c 

Z 



i s - $5 + e - hi 

• • • * 

co = ne - - 3 + + c© 

y 

Equations (A3) and (Ah) are substituted into (Al) and (A2) and then expanded. 
In expanding these equations, some terms will cancel. Terms that involve 
products of three or more angles are dropped from the expansion since they 
represent higher-order effects that are not investigated here. The equations 
for the feathering and flapping moments are then 

M = I (0 + fiB - BC - BC) + (I - I ) [ S2 2 ( 0 - CB) + 0(-B + 20£) - ?B] (A5) 

XX Z y 

M y = i [-g + a(e - d - d) + ?0 + ce] + d x - i 2 )[a 2 B + n(e + b - dWd} 

(a6) 

which define the feathering and flapping moments in the blade frame. The 

inplane moment equation was not retained in this analysis, although the 

inplane degree of freedom is present in the M x and ML. equations. This was 

done to model the forward sweep £ 0 of the blade. Since no inplane motion 

was allowed, the following substitution is made in equations (A5) and (A6): 

C = C=0, C = 0* Also, for later use, these equations are defined in 

terms of their reaction moments M v . . . and M,r. . . , which are given by 

-^inertia v inertia* J 

the negative of the M x and M equations. The equations in the B- frame 

that result from the above substitution are 



hnertia = 'M' 6 ' + ™ ‘ K' ‘ (l Z ‘ V^ 9 ' d 8) “ 


inertia = V® " a{h ~ V )] ' (l x ' + n(9 “ ? o l)] (A8) 


To the flapping equation (a 8) must he added the term due to the effect 
of the hinge offset e. The contribution to the flapping equation of the 
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hinge offset is given by (3/2)l y [e/(R b -e) ]fi 1 2 3 3 (ref . 12) where e is the 
hinge offset and R^, the blade radius. If the effect of the hinge offset is 
included in the flapping equation, the inertia moments in the B- frame are 


M- 


inertia = ' )n20+ Vo M h 


z y o 


(A9) 


^inertia = I ^ 


(A10) 


The above equations will now be transformed to the R- frame (or shaft axis 
system) since the aerodynamic and mechanical moments are developed in this 
reference frame. The transformation between the B and R frames for small 
angles is given by 



(All) 


With the application of the transformation equation (All), the inertia equa- 
tions (A9) and (A10) can be expressed in the R-frame. In expanding the 
equations in the R-frame, the approximation 1 +Cq ~ 1, I x Cq « and drop- 
ping the I x ?o® term in the ^inertia equation (because it was found to be 
small compared to the other terms) gives 


M- 


inertia = ^ ( W h h o Q ' e ~ (l z~\ > ^ e+ ( I x -I kg- ( I « -I ) n B 


+(i -i )c n 2 g- 1 1 z it— 

x y o 2 y o R -e 


(A12) 


^inertia = I b ** 


(A13) 

Equations (A12) and (A13) define the feathering and flapping moments of 
inertia in the R-frame. 


Aerodynamic Moment 

The aerodynamic moment is developed from the normal and tangential 
velocity components at the blade. The assumptions made in developing the 
aerodynamic moments are: 

(1) The aerodynamic center of the blade is at the quarter chord. 

(2) The inflow angle is small. 

(3) Two inflow models were used - a uniform and a nonuniform model. 
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(U) The blade element lift 
coefficient curve is linear. 

(5) The reverse flow region of 
the rotor was taken from 225° to 315° 
azimuth. 

The coordinates used to develop 
the aerodynamic, moments are shown in 
sketch (a). The aerodynamic moment 
about the flapping axis is given by 

r B c 

M = -c l c qr dr (AlU) 

Taero J n 

e s 


In the reverse flow region, the 
limits of integration of equation 
Sketch (a) (AlU) were changed to account for 

the loss of lift in this region. 

The reverse flow region was taken . 
from an azimuth angle of 225° to 315°. 
In this region, the lower limit in the aerodynamic moment equations was set 
equal to the advance ratio y. Although the exact reverse flow region in for- 
ward flight is given by a circular disk of diameter y (for the blade radius 
normalized by R^), this approximation was used since it was more easily 
implemented for computational purposes. To expand equation (AlU), expres- 
sions for the tangential velocity and perpendicular velocity components must 
be obtained. The. dynamic pressure is given by 


1 = 



(A15) 


Since (Up/U^) 2 << 1, equation (A15) becomes 


q - p(u£/2) 

From sketch (a), the angle of attack a of the blade is 


(ai6) 


a = e 


(ait) 


The inflow angle $ is 

<J> = tan - 1 (Up/U^,) (Al8) 

but, for small <f>, a small-angle approximation is 

<f> = u p /u T 

For a symmetrical airfoil, 

c « c = aa 
n Jt 


(A19) 

(A20) 


32 



Substituting for a in equation (A20) yields 


/ U p\ 

; = a 0 - 

n \ V 


(A21 


The tangential and perpendicular velocity components required in equa- 
tion (A21) are given by (ref. 11) 

(r-e) 


Up - ^ 


_ 




(r-e) • 1 

£ + dC cos 1 /) + y sin ifj 

(A22 

• r 1 

8 + — + yB cos tfi - Xj 

(A23 


Neglecting the offset hinge aerodynamically and substituting £ - C G , \ = 0 
yields the tangential and perpendicular velocities: 


T = OT -b( 


P + US cos ip + jj sin tJi 

k D 


■) 


U F = 


3 + TT- 5 B + uB cos </> 
\ ° 


- A ) 


(A2l+ 

(A25 


Nov let 


then 


* = r/f^ 

r = V 

dr ~ R^dx 


(A26 


The aerodynamic moment about the flapping axis can now be obtained by 
expanding equation (AlU). Substituting equations (Al 6 ) and (A2l) into (Aik) 
yields 

“yaero = - ^ f ^ ( BU* - UJjJr dr (A2T 

J e 


PT 


The pitch angle 8 represents the sum of the following angles: 

0 = 0 (x) =* 0 + 9 +x 0 i + 3 tan S 3 

o eye J 


(A28 


where 0 Q is the collective angle of the blade, 0 ^ is the cyclic pitch 

eye 

angle, 0 ! is the twist angle of the blade, and 8 tan S 3 represents the con- 
tribution due to a 63 hinge. Equation (A27) can now be expanded by 
substituting equations (A24), (A25), (A26), and (A28): 
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•V 


aero 


-°p a ^ 2 r B i n 
2 Je /Rb L 


(x+y£ cos ^+y sin ip) 2 (9 +0 + x0 l + 0 l 3 * 1 63 ) 

o o eye 


-(x+yC Q cos \jH-y. sin B+xC^+yS cos *M 0 -* X X cos 4)-X y x sin *)] x dx 

(A29) 

In the above equation, either a uniform or a nonuniform inflow model can 
be used. The inflow X is expressed as 

(A30) 


X = X + X X cos 1 |> + X x sin ifj 
ox y 

where X 0 represents the average inflow and X x and Xy represent a linear 
variation of inflow across the disk. The above equation for X was used to 
integrate the aerodynamic moment equation (My aero ) along the blade span given 

by equation (A29). For the uniform inflow calculation, X x and Ay were set 
to zero. The magnitude of X Q for a given forward flight condition was then 
calculated from the following relationship (ref. 12 ): 

< l/2)C t 

X - yot — — (A31) 

0 S y (y 2 + X 2 ) 1 / 2 


For the nonuniform inflow calculation, X Q , X x , and Ay were estimated from 
measured wind-tunnel data for a particular flight condition. Equation (A29) 
is integrated from the hinge offset location to the tip of the blade (or tip 
loss factor since the equation has been normalized by x = r/R), and making 
use of the definition of the Lock number 


Y 


cpaR^ 


I 

y 


(A32) 


gives, for the aerodynamic moment about the flapping axis, 
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vi a 2 

>c = - -4- 

Taero 2 




cos ip+sin 41 ) 


(5.41 

\ - H)\ 


C o^T “ I|R^ +y ^ cos ^ + ^ 0 sin 40^— | ^ V 2 (Z 0 + Z 0 C0S 2^+sin 2tf>) 

n 4L[C; ..\ 

\ 2 ■ 4J l\t - v 

•4 - k 


|+2y(^ o cos i/i+sin ip) 


cos 2 ip+C Q sin 2 ip 


P-4 

>)(-r - 4, * *■“ • 4 


+2u( ?o cos tp+sin ip)^-£ - - | 00 s 21 ^ sin 2tp^-| _ 3 §tJ 0 1 

+ [(t ‘ TS^ +2 t J (5 o 0os f+sin ip)^-f - 3§^ + 4^f - | cos 2if+5 o sin 2i p\ 

M . 4) 

\ 2 S H/ 


(6 +0 )+ 
o eye 


(“3 " 3lf) +y(? o cos ^ sin *>(t - g|)] X o 
(t - l^ oos nvtcos * sin iP+? o cos 2 <p)^-§ - 

yr - l^j sin *+v(sinVc o cos <p sin ip)^-| - 


To simplify this equation in the main hody of the text, the following 
definitions are made: 


(A33) 


yl ft 2 
K1 * ~|— 


B2 = — - - e l B c _ ®i e 3 
2 2R2 9 B3 ’ 3 ~ 3R 3 


nt b4 e4 ttir - 

Bk = T - TSjT . b 5 = — - 


B 5 e : 


5 5R5 


(A34) 
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GX(tp) = 
C2(*) = 

C3(V>) = 

CM*) = 

C5(*) = 

c6(*) = 

C7(*) = 


5 Bit + y(cos * + 5 sin *)B3 + § W 2 (? 0 + C 0 COS 2* + sin 2*)B2| 


— + 1 (c cos * + sin *)B3 
Q Si o 


B5 + 2y(C o cos ip + sin ip)BU + y 2 ^- - \ cos + ^ 0 sin 2 ^) B ~ 

BU + 2y(& cos ijj + sin vj> )B3 + - \ cos + ^ 0 sin 2 ^) B2 


>(A35) 


B3 + y ( £ cos 4» + sin ip)B2 


BU cos ip + y(cos ip sin ^ + £ q cos 2 iJj)B3 
b 4 sin \f> + y(sin 2 \}> + C cos ip sin ^)B3 


The aerodynamic moment about the flapping axis can now he written in a 
more compact form by use of equations (A3 1 *) and (A35): 

m = -Kl{-C2(i|>) 3 + [ — Cl(i^ ) + C4(ip)tan 63)3 + C3(i|0&i + CU(xJ?)© o 

Taero 

+c 4 (ij))e + C5W\ + c 6 UU + CT(«)x } U 36 : 

eye o x y 

For the reverse flow region, the lower limit of equation (A29) is equal 
to y, which results in a change of the e/R^ terms in equation (A34) to the 
advance ratio y. In the simulation of equation (A36), the appropriate coef- 
ficients were substituted when the retreating blade was in the reverse flow 

region. 

The aerodynamic moment about the flapping axis is given by equation 
(A36). The aerodynamic moment about the feathering axis was found to be very 
small and therefore was not included in the analysis. 


Mechanical Moments 

The inertia moment about the feathering axis of the blade is balanced by 
gyro moments and the flapping moments that couple into the feathering axis 
due to the forward sweep (C 0 ) of the blade. In addition, spring and damping 
moments about the feathering axis must be included. This moment balance is 
established in appendix B where the complete gyro equations are developed. 

For the flapping degree of freedom, the inertial and aerodynamic momen s are 
balanced by the spring restraint of the blade with spring constant k & . The 
mechanical moments about tbe flapping axis can then be written as 
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(A37 


*V Be =h = ye - B p ) 

where 3 and 3p are the flapping and precone angles of the blade. The sum 
of the inertia, aerodynamic, and mechanical moments about the flapping axis 
must sum to zero to obtain a moment balance. The resulting equation for the 
flapping axis is then 

V M Wo- 1 ^ 26 * Z x y ^e- ( i x +i y -i z )!Je-(i z -i y )c o n 2 6-Ki{-c2(v.)l 

+[-Cl(^)+c4(^)tan 5 3 ]3+C3(ip)e 1 +C4(^)e +CU(^)8 +C5(i/j)x +C6(ip)A +C7 (i|j)a } 

o cy coxy 

+k Q { B -8 ) = 0 (A38: 

P P 


This completes the development of the blade equations. The abbreviated nota- 
tion for the aerodynamic terms in equation (A38) was used to simplify the 
equation in the main body of the text. In the remainder of this appendix, 
the flapping equation is developed for any number of blades from which the 
equations for the particular three-b laded rotor are obtained. The flapping 
equation (A38) can be written in general terms for any number of blades as 


Vi-<vyK>“V K sh- 

-Kl|-C2^+(i-l) S2.jB.-Cl[^+(i-l) ^Jg.tan 63 

+C3^+(i-l) S£] 6 1 + Ck[*+(i-l) 2 nJ( 0 o+ 0 cyc ) +C 5 ^ +(i _ l) 

+C 6 [V+(i-l) Sn]x x+ C7[^(i-1) f]x y j+ k B C6i-B p ) -Os i = 1,2,. . .,8 (A3 9) 


where b is the number of blades. The feathering angles in equation (A39) 
are expressed in terms of gyro coordinates 5j and 5 2 . The feathering angle 
0^ is transformed to gyro coordinates and S 2 fry 

6 i = C 2 Js lC os(i-l) ~ + fi 2 sin(i-l) ~~*J ; i = 1,. . . ,b (aUo) 

The derivatives of the feathering angles are obtained from equation (aUo) by 
differentiation. Equation (A39) can now be expressed for a three-bladed 
rotor by varying the index i = 1 , 2 , . . . , b where b = 3 , and by transform- 
ing the feathering angles and their derivatives with equation (A^O). The 
complete flapping equations for the three-bladed rotor expressed in terms of 
the flapping angles and gyro angles are given by 
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i y g 1 -d x+ i y ^-i z )^Bi + | i y eij- W 1 * )nC2 * 1 “( I *" 1 y )5 o o2C241 

I q2 

- — »£ — |-p~ + (c^cos i/j+sin ip)B3^ &-^C 0 B4+ii(oos ifj+r^sin ip)B3 
+ i p 2 (c o +C q cos 2^+sin 2ip)B2j 6i+|b4+2)j ( £ q cos ijj+sin ip)B3+y 2 (|- - ^ cos 2^ 
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+ [B3+yU cos i/>+sin i/>)B2]X + [b4 cos if)+y(cos ^ sin ij;+C cos 2 i|>)B3]X 

o O O -A- 

+ [B^ sin ^+y(sin 2 i|'+C cos sin ip)B3]X |+k (g^-g ) = 0 (A^l) 

o y 1 p 1? 

I ly « 2 M WV^(- I «1 + 4 s *) 
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xL-lz Blt+u[cos(if+120 o )+t sinU+120°)]B3+ | y 2 [? +? cos 2(i|i+120°) 

z lo O 2 OO 

+sin 2(1^+120°) ]B2|6 2 +|BU+2y[ C q cos (iJj+ 120 0 )+sin(if)+120° ) ]B3 

+y 2 [| - | cos 2(4*120° )+e Q sin 2(ip+120° )] B2j g 2 tan <5 3 +|B5+2y [ C q cos ( 4*120° ) 

+sin(ip+120°)]B^+y 2 [| - | cos 2(4*120° )+C Q sin 2 ( 4*120° )] B3J 0j 

+|BlH-2y[c o cos(^+120°)+sin(^+120°)]B3+y 2 [| - ~ cos 2(4i+120° )+C Q sin 2(4*120°)] 

xB2l(e +9 )+{B3+y [ C cos(4*120° )+sin( 4*120° )]B2}X +{BU cos (4*120°) 

J 0 eye o o 


bU+u[cos(i|j+ 120 o )+£ sin(i(i+120 o )]B3+ | y 2 [? +C cos 2(ip+120°) 

) O 2 0 0 


o eye 


+y[ cos ( 4*120° )sin( 4*120° )+C cos 2 (4*120° ) ]B3} X +{B^ sin(i|)+120° )+y[ sin 2 (4*120° ) 
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° y/ p p 


38 



r 7 B 3-c WJ-V° 2 *3 + 1 x y 5^ « 2 fr(yyi z )K 2 (- \ il- 4 « a ) 


~(l -I )c 
z y 


0 ^c 2 (- 1 «i- f « 2 )- (-(% 


+ ^ [^ o cos(4H-2U0 o )+sin{i^+2^0°)]B3 


: $3~ I £ B4+y[cos(^+240 o )+c sin(iH-24o° ) ]B3+ i y 2 [? + £ cos 2(iif+240°) 
V o o d o o 


+sin 2(i|H-240 o )]B2|e 3 + | B4+2y [ C o cos(i|H-240° )+sin(^+240° ) ]B3 

+ u 2 [| - ~ cos 2(iH-240°)+^ o sin 2(ijH-240° )Jb^J 3 3 tan 6 3 +|b5+2u[ C q cos(^+240 o ) 

+sin(i]H-240° ) ]B4+p 2 ^ - ~ cos 2(iH- 240° )+e Q sin 2 ( t/j+240° B3 J 9i 
+ |B4+2y[C o cos (^+240° )+sin(^+240° ) ]B3+ia 2 ^- ~ ~ cos 2(^+240° )+C Q Sin 2 (i|>+240°)J 

xB2 } ( 0 + 9 ) + {B3+y[c cos ( ^+240° ) +sin ( i/j+240° ) ]B2}X +{b4 cos(^+240°) 

f o eye o o 

+p[ cos (^+240° )sin(ip+24o° )+^ o cos 2 (i/j+ 240° ) ]B3>X x + B4{sin(^+240° )+y [sin 2 (?|/+240 o ) 

+C cos(^+240°)sin(^+240°)]B3}X ]+k 0 (3 3 -B ) = 0 (a43) 

u y/ p p 
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APPENDIX B 


DERIVATION OF GYRO EQUATION 


The gyro of the rotor modeled was located at the top of the rotor head 
and was operated at about 10,000 RFM. The gyro and the rotor head are shown 
in sketch (b). This sketch is shown here only for mathematical modeling pur- 
poses and does not necessarily represent the exact mechanical connection of 

the gyro with the rotor. The X s ,Y S , 
Z s coordinates represent the sta- 
tionary coordinates of the rotor. 

Since the body motion of the helicop- 
ter was not modeled, the X S ,Y S ,Z S 
coordinates also serve as the inertial 
coordinates for the development of the 
gyro equations. The X r ,Y r ,Z r coor- 
dinates are the rotating coordinates 
or the shaft axes of the blade while 
the X C ,Y C ,Z C coordinates are the 
coordinates of the gyro at an angle 
i|j 0 from the X r ,Y r ,Z r frame. One 
more coordinate system must be 
defined for the gyro itself since it 
rotates at a much higher speed than 
the rotor. This coordinate system is 
defined by the Xg,Yg,Zg axes system. 
With the definition of these coordinates, the gyro equations and the equations 
for the gyro-blade interaction are derived. 

The rate of change of angular momentum of the gyro is given by Newton's 
second law: 



pV] 

L dt J s " 


M, 


where the angular momentum of the gyro is 

= ^ “tJ/S 


(Bl) 


(B2) 


Since no body motion was allowed (as mentioned in the definition of the coor- 
dinates), the stationary frame X S ,Y S ,Z S serves as the inertial frame in the 
development of the gyro equations. 


The rate of change of angular momentum of the gyro with respect to the 
frame is expressed as 


d0) G/S 
dt ( 


+ “c/s * x e 


Vs = m g 


(B3) 


UO 



To expand equation (B3) , the angular velocities "between reference frames 
must "be established. The angular velocity of the gyro with respect to 
inertial space is 

Vs = Vc + “c/R + “R/S (Bl+) 

The angular velocity of the C frame with respect to inertial space is 

“c/S = “c/R + “R/S (B5) 

The angular velocities in equations (b 4) and (B5) are 


J C/R 


= Csi 5S2 5S3 1 




0 
0 

LMJ 

Ri 

62 

LDJ 

’0 

0 

lP-I 


(b6) 


(BY) 


(B8) 


The transformation between the C and G frame, considering only small 
angles, is 


C 1 


1 

0 

-5 2 

r-i 

Sl 

C2 

- 

0 

1 

6 1 

g2 

C 3 


6 2 

-Si 

1 

S3 


(B9) 


Equation (b 4) can be expressed in the C frame by use of equations (B6), 
(B7), and (B8) and the transformation equation (B9): 


hi 


G/S 


~ tci.c 2 .c 3 ] 


6 1 — £2 6 2 


(BIO) 


where oiq. = o> + ft 
The angular rate 
differentiation: 


and 

W G/S 


• * 

terms such as 8 i 5 2 and 5 i 6 2 are considered negligible, 
can be obtained from equation (BIO) by 


hi 



“ [ C 1» C 2> C 3^ 


n 

•« • 

— ft 62 


^2 + ft<$i 
0 


_l 


(Bll) 


Equation (B5) can also be expressed in the C frame by use of equations 
(B7) and (B 8 ) and transformation equation (B9): 


^C/S = 


6 j — ft 6 2 
^2 + ftS]_ 

n 


(B 12 ) 


Since the gyro is modeled as a disk, the diametrical inertias are equal, 
•defined as I G = Ij = I 2 - The inertia through the center of rotation of the 
disk is I 3 . °The inertia distribution of the gyro is then 



w — 


— - 


X G 0 0 


Cl 

[c\ ,C2,C 3 ] 

0 I G 0 


C2 


0 0 I 3 


C3 

mm 


(B13) 


Substituting the above equations for the angular velocities, rates, and 
inertia distribution into equation (B3) gives the following equations for the 


gyro; 

I^6i + ( X 3^2 o> g - I^ft 2 ) 5 1 + (i 3^ - 2 fiI G )62 = (BlU) 

I g '« 2 + (X 3 Qo> g - I g S2 2 )S 2 - (l3» 0 - 2BI g )&i = (B15) 

For a high-speed gyro, will be some multiple value of the rotor 
speed ft. Hence let 

U G = (Bl6) 

where kp is some multiple of the rotor, speed. Substituting equation (Bl6) 
into (Blfi) and (B15) gives for the gyro equations: 

I G «1 + (I 3 n 2 kj, - I G n 2 )«i + (13^ - 2ni G )«2 = (217) 

X g «2 + (l 3 n2 kfi - i g « 2 )«z - (13^2 - 2ni 0 )ix = MQy (B18) 


In the following section, the external moments Mh and Mq of equations 
(B17) and (Bl 8 ) are developed. ^ 


U 2 



Z», , Z 


c 



'Vs = ^2 - C s e , 


First, the pilot 1 s svashplate 
roll and pitch moment inputs (vj 
and v 2 ) and the svashplate damping 
and spring terms are developed and, 
second, the excitation of the gyro 
by the blade is derived. The sum 
of the above moments in rotating 
gyro coordinates constitutes the 
external moments Mg. and Mg . A 

y 

sketch of the gyro disk ■with the 
springs , dampers , and moments 
inputs is presented in sketch (c). 
These moments were measured in the 
stationary axis system (X S ,Y S ,Z S ) 
and must be transferred to the 
rotating gyro axis system (X^Y , 

z 0 >. 


The moments in the stationary 
axis system due to pilot, input and 
the spring and damper are 


- K 4> 
s sp 

(B1 9 ) 

- k e 

s sp 

(B20) 


The above equations must now be transformed to rotating coordinates since all 
equations in this development are expressed in the rotating coordinate system. 
The transformation between the stationary and the rotating coordinates is 
given by 



sin(i|H-i|0 

- 

cos(ijH*\p ) 

Cl 

-cos(i/>+^ ) 

sin(t|>+ij; ) 

o_ 

c 2 


(B21) 


or, between gyro angles and svashplate angles, by 


4 > I r sin(i]j+iii ) 
sp _ r r o 

<f> -cos(iH-^ ) 

L spj L o 


cos(iJ;+iJ> ) 
s in ( ) 



(B22) 


Equations (B19) and (B20) can be transformed to rotating coordinates by apply- 
ing transformations (B2l) and (B22). In addition to the pilot input and the 
spring and damper moments, the excitation of the gyro by the blade 

exc^c 

and Mexcyc must "be added to give the complete set of applied external 
moments for equations (BIT) and (Bl8): 


■U3 


(B23) 


M Gx = ^sinfip + H> Q ) - v 2 cos(ifj + ^) + ftC g 6 2 - C g 6i - K g 6 i + M exc X c 
KU v x cos(^ + ^ Q ) + v 2 sin(i/j + rp Q ) - fiC^ - C g 6 2 - K g 5 2 + ^exc yc 

(B2U) 

The gyro equations can be written from equations (BIT), (Bl8), (B23) , and 
(B2^) as 

I G 6l+(l3ft 2 ^-lQ^ 2 ) ^ 2 ”‘^C s ^ 2 + ^ s ^l + ^- s ^ i^excxc 

= v 1 sin(^+^ ci )-v 2 cos (^+ip o ) (B25) 

I G 5 2 +(l3«\-I^ 2 )52-(l3^V 2fiI G )Sl+fiC s 6l+C s^ 2+K s 62_Mexc yc 

= v^cos (^+V J 0 ) +v 2 s i n ^ +l J J, 0 ) (B2 6) 

The only remaining terms to be expanded is the excitation of the gyro by 

blade M exc and V^xc . These are the moments in gyro rotating coordinates 
xc yn 

(X ,Y ) which are transmitted through mechanical linkages from the feathering 
axis of the blade. The moment about the feathering axis for one blade in the 
rotating blade coordinate (X r ,Y r ,Z r ) system is given by (ref. 13) 


M 

exc 




mech 


C o ] V 


inertia 


+ M. 


x inertia 


- V 


V 


(B2T) 


The moments ^y^^ and ^inertia are ' the mechanlcal inertia moments 

that couple into the feathering axis because of the forward sweep of the 

blade. The inertia moment about the feathering axis is given by 


M x . + . . 

inertia 

The kgQ and Cq© terms represent the spring and damping effects about the 
feathering axis. Equation (B27) is expanded by substituting for the inertia 
and mechanical moments from equations (A12), (A13), and (A3T ) . In the expan- 
sion, the approximation is made that 1+C 2 - 1 and I X C Q is small compared to 
other terms in the equation: 


m = -i e-c a e-(i -i )^ 2 e-k e-(i +i -i )ri+(i -i k fl 2 m fi ( 6-6 ) 
exc x 0 2 y 0 x y z zyo o0 p 


(B28) 


The moment equation (B28) is expressed in the rotating system of the blade 
(Xj.jYp) and is now transferred to the gyro axis system (X C ,Y C ). The accom- 
panying sketch (d) shows the three-bladed rotor treated in this development 
with the rotating frames of the blade (X r ,Y r ) and the gyro (X C ,Y C ). 

The moment M exc about the feathering axis of each blade (sketch (d)) 
will be resolved along the (X r ,Y r ) axis and then transformed to the (X C ,Y C ) 



For word 

t 



Sketch (d) 


axis through a mechanical linkage constant C 2 - The excitation 1 of the gyro 
by the blade is then given by 


= c 2 m (ai,e 1 )-.~M (62,62) -i-M (63,63) 

excxc exc 19 1 2 exc Z9 z 2 exc 39 3 J 

Mgv. = C2|"^ M (62,82) - ^ M (63,63) 
excyo ^ 2 exc z z 2 exc 59 3 I 

Substituting for each moment in equations (B29) and (B30) from (B28) with the 
appropriate feathering and flapping angles for each of the three blades gives 
the following moment equations: 



- k e( 0 i- I ® 2 - | 1 82- | 63) 

+( V I y k o S52 (^- I ®2- I ® 3 ) + ^o k 6 [( 8l“6p)- | (6 2 -6 p )- § (B S -B p )]} 


(B31) 


(B29) 

(B30) 
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= 41 c 2 {-i (S 2 - 03 )-o n (e 2 -e 3 )-(i -i )n 2 (e 2 -e 3 )-k (e 2 -e 3 ) 


' yc 2 


-(l x +I y -I z )S2(B 2 -g 3 ) + (l z; -I y )C o n 2 (P2-B3) + C 0 k B [(S 2 -e p )-(B3-S p )]} 


The feathering angle 0^ (i = 1,2,3) and its derivatives are now written in 
terms of the gyro coordinates 5^ (i = 1,2). The relationship "between 0^ 
and 6^ and the derivatives can he obtained from equation (A^O). Applying 
this transformation to equations (B3l) and (B32) gives the excitation of the 
gyro by the blades in gyro coordinates as: 

M exc x0 = - fctl x V f c Cl*!- | c|(l J 1 - fv&l 


2 z V 


-C 2 (l x +I y -I z )S2($ 1 - | B 2 - \ B 3 )+C 2 (X z -I y )? o fi 2 (g 1 - § B 2 - \ g 3 ) 
+ C 2 k 3 [(B 1 -B p )- | (B 2 -B p )- § (B 3 -B )] 


Mexc yc = ' | C|l x « 2 - |c 0 cl5 2 . | c| ( I z _I y ) S T 2 ^5 2 - fkgC^ 

- 4 C 2( I x + Vb )n( ' 62 " S3)+ ^ C 2 Cl z -I y )C o n 2 (B 2 -B3) 

+ C 2 k 6 [(6 2 -B )-(6 3 -B )] (B3l*: 

Equations (B33) and (B34) are combined with equations (B25) and (g2 6 ) to pro- 
duce the complete gyro equations for a three-bladed rotor: 

(v f C 2 I x )«l + | o e c 2 Bi+c s 6 1 +(i 3 Q 2 k B -i G n 2 )6 1 + |c!(i 2 -x y ) n 2a 1+ |k 0 C 2 « 1+ K s Sj 
+ ( 1 3 sik R -2ni G ) a 2 -nc s 5 2 +c 2 ( i x +i y -i z ) n ^B i - \ B 2 - f s 3 )-c 2 (i z -i y k o n 2 
x (bi- t I B 3 )-C 2 C o k g [(B 1 -B p )- \ (B 2 -B p ). | (B 3 -a p )] 

= visinC^+ip )-v 2 cos(’B+' 3 | o ) 


(B35) 



c e c|6 2 +C g 6 2 + (l3« 2 V« 2 I G )«2 + |cl(l z -I y )« 2 *2 + |k 9 C2«2««2 

-(l 3 fik H -2S2I G )S 1 +nC s « 1 + ^ C 2 (l x +I y -I z )S2(ES 2 4 3 )- ^ c 2 (i z -i y )t; o n 2 (B 2 -B3) 

jy C 2 q o kg[ (B 2 -B j) )-(B 3 -B )] = v 1 cos(i{i+iii o )+v 2 sin(i|)+')» 0 ) (B36) 


(V I C2l x)^ + I 
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APPENDIX C 


hub moment calculations from measured data 


This appendix develops the transformations necessary to change the mea- 
sured shaft "bending moments and load cell moments to moments at the huh. The 

shaft "bending moments and load cell 
moments were measured during wind- 
tunnel testing at some distance "below 
the hub; the moments obtained from 
the math model simulation were cal- 
culated at the hub. Therefore, it 
was necessary to transfer the mea- 
sured moments to the hub in order to 
have a common point of comparison 
between measurement and simulation. 
The coordinate system used to cal- 
culate the hub moments is shown in 
sketch (e). 


'LC, 







Sketch (e) 


4 = 


M- 


LC 


x 




The load cell moments were mea- 
sured in stationary coordinates a 
distance h^c "below the hub, while 

the shaft moments were measured in 
rotating coordinates a distance 

below the hub. In calculating the 
moments at the hub in rotating coor- 
dinates, the following notation and 
positive direction of moments are 
adopted: 


(Cl) 


rotating 

load cell roll moment (positive left roll) 


load cell pitch moment (positive nose-up) 


S- 






rotating 

shaft bending moment (positive left roll) 


shaft bending moment (positive nose-up) 


(C2) 


The desired hub moment in rotating coordinates is defined as 
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%; 

rotating 

hub roll moment (positive left roll) 



hub pitch moment (positive nose-up) 


The force acting through the huh is defined as 


" F B X " 

rotating 

hub force along 



hub force along Y^ 


(C3) 


(CU) 


The following moment balance can be written from sketch (e): 

"2c - ^ + e lc x *b 

# = + h g x F« 


(C5) 

(C6) 


Equations (C5) and (c6) can be solved simultaneously for the. moments at 

the hub from the known measurements of 52* . M R , hL_. and h . The hub moments 
. .. . „ , XC S s a LC* s 

are then given by 


Mi 


h LC z Mx r - h s, M LC. 


h LC z - h Sz 


% = 


h LC z *V r - h s z MLC y 
h LC, " h s„ 


(C7) 

(C8) 


The load cell moments expressed in equations ( CT ) and (C8) are in the 
rotating coordinate system. Since the. measured load cell moments were 
recorded in stationary coordinates, the transformation between the rotating 
and stationary coordinate system must be obtained. This transformation, as 
indicated in sketch (e), is given as 




— ^ 


— — 

*1 


cos t/i sin ip 0 


si 


Si 

-sin ip cos ip 0 


S2 



0 0 1 


S3 

— — 


— _ 




(C9) 


Applying the transformation given by equation (C9) to the and M^q 

terms in equations (C7) and (c8) produces the final form of the hub moments 
in terms of the measured moments: 


y 



h LC g M x r h s 2 ( M LC x cos + M LC y sin ^ 
Rg x ' h L C z - h s z 

h]X.My r - h Sz (-M L g x sxn » + M LCy cos ») 
% " h LC z - h s z 


(CIO) 


(Cll) 



Sketch (f) 


In the mathematical model simu- 
lation of the rotor, the huh moments 
were calculated from the blade flap- 
ping angles and blade spring con- 
stant. Sketch (f) shows the 
coordinates used to obtain the hub 
moments in the rotating coordinate 
system. From the geometric location 
of the blades for a three-bladed 
rotor, the hub moments are expressed 
as 

M Sx = * 3 (t - 4 03 ) (C12) 



(C13) 


The calculated hub moments given by equations (C12) and (C13) were then 
compared with the measured hub moments expressed by equations (CIO) and (Cll). 
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□ 

n(iJi+120 o )+? cos 2 (i{-+120 o )]B3}X 


C cosd+120°)sind+120°)]B3}A +7^0 
0 y r y p 


B3+u 2 ^ - cos 2('Jj+240°)+r > ^sin 2 (iJj+24o°)J B 2U 
)]B4+y 2 j4- 4 COS 2( 1(1+240° )+C o 3 in 2(i|»+240° )J Bsj 9 j. 
]B2}X 

O 

n(?(.'+240°) + C cos 2 (iJ>+240 o )]B3}X 


1 1 

.+ 7 y^ 2 (B4 sin(i(i+240 a )+v[sin 2 {4i+240 o )+c cos(ij/+24o 0 )sin{<ii+24o 0 ) ]B3lX + r~ B 

y I y P 
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TABLE 2«- ROTOR AND GYRO PARAMETERS 


Parameters 

Nominal values 

with allowed percentage variation 

Rotor: 


Number of blades (b) 

3 

Blade inertia (for each of 


b blades ) 


Feathering inertia (l x ) 

0.54 kg-m 2 (0.4 slug- ft 2 ) 

Flapping inertia (l^) 

363.4 kg-m 2 (268 slug- ft 2 ) ± 1 percent 

Inplane inertia (l 2 l 

363.9 kg-m 2 (268.4 slug- ft 2 ) ± 1 percent 

Blade twist (Qj) 

-9.43° ± 5 percent 

Collective angle at root ( 0 Q ) 

8 . 57 ° ± 7 percent 

Sweep, angle ( Cg) 

1 . 5 ° ± 20 percent 

Precone angle (6 ) 
P 

2 . 25 ° ± 5 percent 

Tip loss factor ( B ^ ) 

0.97 ±2 percent 

Blade radius (R^) 

5.03 m (16.5 ft) 

Blade chord width (c) 

0.357 m (1.17 ft) 

Airfoil type 

NACA 63 2 015 

Solidity 

0.0675 

Blade area 

5.36 m 2 (5T.T ft 2 ) 

Disk area 

T9.U6 m 2 (855.3 ft 2 ) 

Mast angle forward tilt 

6 ° 

Rotor RPM ( 100 percent ) 

355.6 

Hinge offset (e/R-b) 

0.11 ± 30 percent 

cpaR^ 


Lock number, y = — = — 

V 

4.57 ±10 percent 

Spring constant (p from fig* 4) 

±30 percent 

k = I Q 2 (p 2 -1- -^-- e - f- tan 


3 y l R-e 8 y 


Gyro: 


Gyro inertia 


Polar moment of inertia (l 3 ) 

0.4l kg-m 2 (0.3 slug- ft 2 ) ± 2 percent 

Diametric moment of inertia (ip) 

0.21 kg-m 2 (0.15 slug-ft 2 ) + 2 percent 

Gyro pitch link angle (^ Q ) 

60 ° ± 2 percent 

Mechanical advantage constant (C 2 ) 

0.87 ± 3 percent 

Gyro RPM (100 percent) 

10,000 


preceding page blank not 


film® 
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TABLE 3.- PARAMETER ESTIMATES FOR STEADY-STATE LOCKED GYRO OPERATION 


Condition 

Parameter 

Value 

Nominal 

Estimated 


6 o 

8.57° 

7.97 


B £ 

0.97 

0.95 


C 2 

0.87 

0.85 


Iy 

363.1+ kg-m 2 

359.7 kg-m 2 

V = 1+1.2 m/sec 

y 

(268 slug-ft 2 ) 

( 265.3 slug-ft 2 ) 

(80 knots ) 

i 2 

363.9 kg-m 2 

360.2 kg-m 2 

y = 0.399 


( 268 . 1 + slug-ft 2 ) 

( 265.7 slug- ft 2 ) 

rpm = 55.5 percent 

Y 

lw67 

5.13 



ll+l s 092 Nm/rad 

129,31+5 Nm/rad 


3 

(I0l+,061+ ft- lb /rad) 

(95,1+00 ft- lb/rad) 



1.5° 

1.6° 



2.25° 

2.36° 


x 0 

-0.012 

-0.012 (fixed) 


e/Bb 

0.11 

0.11 


e o 

8.57° 

7.97° 


B £ 

0.97 

0.95 


c 2 

0.87 

0.85 


I 

363 .H kg-m 2 

359.7 kg-m 2 

V = 31 „ 3 m/sec 

y 

(268 slug-ft 2 ) 

( 265.3 slug-ft 2 ) 

(6l knots) 

i z 

363.9 kg-m 2 

360.2 kg-m 2 

y = 0.1+02 


( 268 . 1 + slug-ft 2 ) 

(265.7 slug-ft 2 ) 

rpm = 1+1.6 percent 

Y 

1+.69 

5.16 


k R 

125,211+ Nm/rad 

87,651 Nm/rad 


3 

(92,353 ft-lb/rad) 

(64,648 ft-lb/rad) 



1.5° 

1 . 6 ° 



2 a 25° 

2 . 1 I + 0 


X o 

0.005 

0.005 (fixed) 


e/R b 

0.11 

0.11 
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TABLE 4 . - ' PARAMETER ESTIMATES FOR TRANSIENT FREE AND LOCKED GYRO OPERATION' 
v 55 3 6 m/sec (TO knots), y = 0.355, 54.2 percent npm 



Condition 


Values 


of gyro 

Parameter 

Nominal 

Estimated 



e o 

8.57° 

9 . 17 ° 



01 

-9.43° 

-10.93° 



% 

0.97 

0.96 


Free 

Y 

>*.57 

U.O 



h 

138,752 Nm/rad 

97,127 Nm/rad 



(102,338 -ft-lb/rad) 

(71,637 ft- lb /rad) 



bo 

1.5° 

1.8° 



e o 

8.57° 

9.17° 



B £ 

0.97 

0.98 



I 

363. 4 kg-m 2 

366.5 kg-m 2 




(268 slug-ft 2 ) 

( 270.3 slug-ft 2 ) 


Free 

Y 

4.57 

lt.0 




138,752 Nm/rad 

97,127 Nm/rad 

Parameter 



(102,338 ft-lb/rad) 

(71,637 ft- lb /rad) 


Co 

1.5° 

1.8° 

cases 

studied 


a *b 

0.11 

0.15 



®o 

8.57° 

9.17° 



B it 

0.97 

0.98 



C 2 

0.87 

0.85 



I 

363.4 kg-m 2 

359.7 kg-m 2 


Locked 


(268 slug-ft 2 ) 

( 265.3 slug-ft 2 ) 


Y 

b. 57 

5.02 



H 

138,752 Nm/rad 

97,795 Nm/rad 



(102,338 ft -lb /rad) 

(72,130 ft-lb/rad) 



V 

0.355 

0.34-8 



Co 

1.5° 

1.8° 



3 

p 

2.25° 

2.1V 3 



e/R b 

0.11 

0.15 
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TABLE 5.~ SENSITIVITIES OF MODEL PARAMETERS 
OBTAINED FROM IDENTIFICATION PROGRAM 


Parameter 

Sensitivity 

constant 


4.0 H 

e l 

3.2 

% 

2.7 

c 2 

1.1 

* ( VV 

.78 


.75 

V 

.72 


.57 

I G 

.57 


.52 

% 

.31 

x o 

.15 

c s 

.13 

c 0 

.07 

I x 

.03 

£ 

.02 
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